moagem controlada até ao calibre desejado constituirá um incremento divisível sem perda 
de representatividade que, com relativa facilidade, permitirá satisfazer aquelas condições. 


— Mas, se os incrementos houverem que ser colhidos em galerias, etc., terão de ter peso mínimo 
para representarem a granulometria que exibirão e podem, por isso, ser relativamente volu- 
mosoz. Mas, mesmo assim, se o processo de tratamento a ensaiar exigir, a priori, moagem 
abaixo de certo calibre fino, poderão ser controladamente moídos e divididos de forma a cons- 
tituirem convenientemente os pequenos incrementos desejados. 


(ii) Todavia, pcderá dar-se o caso, aliás frequente, do próprio processo de tratamento em investiga- 
ção desejar esclarecer o comportamento do minério em calibres relativamente graúdos. Nesse 
caso, há que abandonar, quase sempre, o recurso aos testemunhos das sondagens e utilizar como 
locais de amostragem as galerias, etc. e ndo será possivel constituir incrementos de pequeno 
porte que sejam representativos dos inter valos AZ, arbitrariamente escolhidos. 


Os incrementos terão de ter peso que se não poderá fixar a priori, sem se considerar 
as exigências da teoria de Gy. 


2.º -— A acessibilidade à jazida é reduzida. Por exemplo, os testemunhos das sondagens não existem, ou 
não são de confiança, ou já se alteraram, não existe ossatura de preparação suficientemente extensa 
ou, a que existe é tal, que só em certos locais se podem colher amostras; ou, ainda, há que abrir gale- 
rias, travessas ou chaminés, que terão de ser de reduzida extensão, porque o tempo é escasso ou, 
também, não é possivel, por razões econômicas ou de temporização, realizar conveniente campanha 
de sondagens adicional. 


Nesse caso, a escolha dos locais de amostragem estará confinada apenas a algumas possibili- 
dades, 

Suponhamos então que se dispõe de um ponto acessível que se localiza num bloco ou painel V. 
se ai colhermos um incremento, a variância de extensão, isto é, o estimador da variância do erro de 
atribuir à amostra o teor krigado do bloco ou painel V é, como se sabe, 


a , * 
1 j Es 
| |- dv ta | ds | 7 (u—y) dy = 04, + Y(V,V) 
= r Y 


a v a 


Poderá, por exemplo, dizer-se que haverá aproximadamente a probabilidade de 95º para que o teor da 
amostra se afaste do teor estimado do painel mais de2y' elo 


Se a situação for um pouco mais favorável, isto é, se, em dado painel ou bloco, se puderem colher 
vários incrementos (sejam N), o estimador da variância de extensão ao teor médio desses incrementos 
do teor estimado para o painel, é, obviamente, 


a * E) E es INT ATA l à 


VN 


1 wu 


Poderá, portanto, procurar-se, dentre os locais acessíveis, aqueles que revelem maior probabilidade de 
fornecer incrementos com certos teores médios, previamente escolhidos para preencher o histograma 
a reproduzir. 


Mas, seja como for, colhidos os incrementos, eles exibirão individualmente certos teores médios 
que, por sua vez, representarão nova população de amostras que há que comparar com a de referência, 
atrás mencionada. 


Quer dizer, nestes casos de reduzida acessibilidade da jazida, não há, em geral, forma de repro- 
duzir o histograma de referência — obtem-se uma população de incrementos que, eles próprios, exibem 
certo histograma, em geral, distinto do de referência. 


O primeiro passo a dar é o da comparação dos dois histogramas (o dos incrementos e o de refe- 
rência) para verificar se dizem respeito à mesma população. São estes clássicos. Com muita frequência, 
estes testes são negativos, mesmo tomando o cuidado de relacionar o histograma de referência a suporte 
da ordem de grandeza do dos incrementos ou vice-versa, 
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Assim sendo, um dos caminhos possíveis num caso destes, é construir amostra compósita, com 
os incrementos Z, que satisfaça às condições: 


b'(Z — m) mto () a br = À 
(D(Z)—D(Z,)P— O 4» = 0 
1 E [1,5 NJ 


Se for possível satisfazer I com cs pesos dos incrementos e os respectivos teores disponíveis, 
ter-se-á amostra única compósita representativa, quanto ao atributo em questão, no que diz respeito 
à média e com dispersão de teores que se afasta o menos possível da da distribuição de referência. 
Quando houver interesse simultâneo em vários atributos X, Y, Z, etc., poder-se-á compor, com os incre- 
mentos, amostra que satisfaça a condição muito mais forte: 
[b'|X,-E(X)|— 0 
mil = [9:*)] = 
[E (o Dib l 
i 
| (Te = Ea a 
' | b'[Z,-E(Z)|-—o0 
| 


Er [DX )-D (XD PA+s[DOYy-DODP +He—0 de [Ema 
| r4+s+t+.=lirs,..>o 


Condição que só com muitos locais acessíveis será possível satisfazer, como é óbvio, 

É claro que, nestes casos, haverá uma amostra única, cujo pezo, muitas vezes, não será diminuto, se 
se atender à granulometria com que os incrementos são obtidos. Dependerá do tipo de teste de tratamento 
a realizar. 

Todavia, nem sempre sistemas como I e I' serão trabalháveis. Poderá então comparar-se a popu- 
lação de incrementos com a de referência em termos de média e de variância. 

Esta via pode ter mais interesse quando há, simultaneamente, vários atributos a estudar. 


Exemplo (Figs. 1, 2 e 3) 


O exemplo diz respeito a amostragem realizada em extensa bacia carbonifera Brasileira relativamente 
regular, onde se havia definido, por numerosas sondagens, duas camadas sobrepostas de carvão e onde a 
acessibilidade era muito reduzida. 

A bacia havia sido, arbitrariamente, dividida em quatro áreas denominadas I, II, HI e IV. E as cama- 
das classificadas de «superior(s)» e «imferior(i)». 

Foi possível colher incrementos nessas duas camadas apenas na área II, que constituiram dua: popu- 
lações denominadas por A(II,s) e A(ILi). 

A comparação entre os histogramas de A(Ils) e A(ILi) com os histogramas de referência de qual- 
quer das áreas (histogramas deduzidos dos testemunhos das sondagens) mostrava que os teores de cinzas, 
matérias voláteis e enxofre não pretenciam às mesmas distribuições embora os suportes fossem aproximada- 
mente os mesmos. 

Optou-se então em estudar as equivalências em média e em variância das populações para, de certo 
modo, classificar a representatividade das duas amostras colhidas, As figuras referidas mostram os resul- 
tados obtidos, os quais permitem, sem dúvida, algumas conclusões pertinentes quanto a representatividade. 


3.2. Segundo objectivo: Afinação (optimização) do processo de tratamento 
Construção de curvas de lavabilidade semi-industriais 


Na realidade, para cumprir este objectivo, o suporte «bloco tecnológico» e respectivo histograma pode 
ter reduzido interesse. Uma vez que o processo de tratamento a adoptar tenha sido escolhido, o objectivo da 
amostragem é então o seguinte: «Pretende-se regular (e optimizar) o processo para uma alimentação que 
deve simular a futura alimentação do concentrador. Com frequência, a construção da curva de lavabilidade 
semi-industrial poderá estar incluída nos objectivos desta investigação final.» 

Assim sendo, as populações dos incrementos devem ser representativas, cada uma, de cada tipo de 
minério que haja sido definido na jazida, o que não traz novidade relativamente ao que se referiu, mais 
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como a maioria dos processos se realiza efectivamente sobre produtos de calibre relativamente fino, o que 
interessa é considerar os histogramas «pontuais», 

Há, portanto, que tentar realizar krigagem para suportes «mínimos» e construir os respectivos histo- 
gramas que serão mais dispersos do que os das carotes. 

Interessa então simular amostragens referidas a estes histogramas, usando, por exemplo a metodologia 
atrás mencionada, 

Há aqui, todavia, outra vantagem adicional, é que se lida normalmente com grandes amostras, pois 
se realizam testes piloto ou semi-industriais, o que implica maior acessibilidade à jazida. 


MATRIZ DE EQUIVALÊNCIAS GRAFO DE EQUIVALÊNCIAS 


S(IlI,s) SUIL,i) S(III,s) S(ILI,i) S(IV,8) A(lI,s) A(IL,i) 


S (11,5) 


S(IL,s) 1 — 1 — 
S(II,i) — 1 -— 
S(Ill,s) 1 — 1 
S(MIji)j| — 1 Ui) 
S(IV,s) LI) — O 
A(II,s) 1 1 
A (II) — 1 — 1 — 


N. B. 1— equivalente; O — não equivalente. 


o (II, 5) 


| 
| 


E E E o. = + 


S (II, i) 4 All, i) 


equivalente 


-———— não equivalente 


S (Hi) 


Fig. à 


Equivalência em média quanto a cinzas 


MATRIZ DE EQUIVALÊNCIAS GRAFO DE EQUIVALÊNCIAS 


S(LL8) SLi) SUL s) SOIL) SUV,8) AML) A(ILÃ) 


S (1,5) 


S(II,s) 1 — 1 — O 1 -— 
S(II,i) — 1 — 1 O O 1 
S(II,s) 1 — 1 ) 0 1 — 
s(WIi)| — 1 0 1 E A (1,5) 
S(IV,s) O — O — 1 0 — 
A (II) 1 — É - 0 1 1 
A (II) — 1 — E — 1 1 
N. B. 1 — equivalente; | O — não equivalente. 

S (Ui) Adu, i) 

equivalente 
.—-—-—-- Não equivalente 
S (11,1) 
Fig. E 


Equivalência em variância quanto a cinzas 
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GRAFO DE EQUIVALÊNCIAS 


S (11,5) 


Fig. 3 


— equivalente 


====== não equivalente 


Equivalência em média e variância quanto a cinzas 


4. LISTA DE SÍMBOLOS NÃO DEFINIDOS NO TEXTO 


Z atributo, característica, variável regionalizada 
Z.. valor médio do atributo no volume Y (bloco, painel, bloco tecnológico) 


Ar estimador de krigagem de 4, 


E( .) esperança 
D(.) variância de dispersão 
A A intervalo de atributo 


AZ valor médio arbitrado para o intervalo 4 Z 


Y(.) variograma 
y(.,) variograma médio 
m valor médio = E (.) 


(Continuação da pdg.º 208) 

A British Steel Corporation que, 
consoante a política vigente, so- 
freu os efeitos de uma naciona- 
lização (1949), desnacionalização 
(1958) e mova nacionalização 
(1967), procedeu a uma redução 
suplementar de efectivos (38 mil). 
Por seu turno, a ARBED, impor- 
tante empresa do Sarre, que, com 
a tomada de controlo da Neunkir- 
chen na RFA e da Rodange-Athus 
no Luxemburgo, pôs ponto final na 
reestruturação da siderurgia regio- 
nal, deverá despedir 8200 efectivos, 
até 1982, para poder atingir o 
ritmo de rendibilidade dos japone- 
ses. 

No caso de pretender adaptar o 
seu aparelho de produção aos no- 
vos dados do mercado do aço, ou 
seja, ao abrandamento, que tudo 
leva a crer seja duradouro, do 
ritmo de crescimento do consumo 
(consequência retardada dos efei- 
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variância de krigagem de Z,, 


tos da crise sobre os clientes das 
siderurgias), e à entrada de novos 
paises produtores no mercado (a 
produção de aço brasileira registou 


Quatro anos de «crescimento negativo» 


Paises 

Grã-Bretanha. 

Itália 

França 

Rep. Fed. da Alemanha 
Bélgica 

CEE 

USA 

Japão 


| 975977. 
2 1975-1976 


um aumento de 229% entre 1960 
e 1974), a siderurgia europeia terá 
ainda de solucionar sérios proble- 
mas, 


Evolução de 1974 a 1977 em % 


Taxas de Produção Produti- 
utilização da vidade 
capacidade 
— 12 — 10 — 1 
— 15 — (0,2 — 8 
— 22 — 22 —17 
— 30 — 36 — 26 
— So — dd mim A 
—24 78 —16 
+ d! —17 — 1 
—15? —14 —10 


[Extraido do «Boletim Mensal BFN», 
(12) Dezembro 19781]. 


(Continua na pdg.º 260) 
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Introdução 


RESUMO 


Os principios da simulação estocdástica são dis- 
cutidos na sua generalidade e o Método de Monte 
Carlo é descrito brevemente. 

Aplicações da Simulação Estocástica nos dife- 
rentes domínios das Ciências da Terra são referidas 
no texto e bibliografia existente no CVRMUL é indi- 
cada para cada caso especifico. 

Os domínios considerados são a Geologia e Pros- 
pecção, Exploração e Tratamento de Minérios e fina- 
liza-se com a Análise económica global de um pro- 
jecto mineiro. 
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à Simulação 


HENRIQUE GARCIA PEREIRA 
Instituto Superior Técnico 


Universidade Técnica de Lisboa 


ABSTRACT 


The principles of stochastic simulation are briefly 
discussed and Monte Carlo method is described. 

Applications of stochastic simulation in some 
fields of the Earth Sciences are referred and referen- 
ces are given for each specific case. 

The main fields treated are geology and explo- 
ration, mine planning and ore dressing. The global 
economic study of a mining venture is also introduced 
in a case story. 


1. GENERALIDADES 


A simulação é uma técnica cujo objectivo é re- 
produzir o comportamento de um sistema complexo 
através de um modelo, 

Esta técnica é em geral aplicável quando a expe- 
rimentação exaustiva sobre o sistema é difícil ou 
proibitiva do ponto de vista económico. 

Utiliza-se ainda a simulação em sistemas em que 
a base teórica que descreve o seu comportamento no 
tempo (por exemplo, sistema de equações diferen- 
ciais) é completamente desconhecida ou em que cer- 
tas operações unitárias são conhecidas, mas as inter- 
dependências entre estas são complexas ou não 
formuláveis de um modo coerente, 

Neste tipo de problemas, é usual recorrer à simu- 
lação do sistema através da construção de um modelo 
que descreva na sua complexidade e detalhe a reali- 
dade que se pretende reproduzir (e não necessa- 
riamente na sua «essência», como se pretende no 
«método científico»), 

Na base de toda a simulação está pois a formu- 
lação de um modelo, e esta vale na medida da vali- 
dade do modelo. Sem nos determos agora sobre a 
construção dos modelos (admitimos que estes repro- 
duzem o sistema numa determinada gama de condi- 
ções), vamos apenas classificá-los em FÍSICOS (me- 
cânicos e analógicos [1]) e MATEMÁTICOS. 
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Trataremos apenas da simulação executada so- 
bre modelos matemáticos. Estes são, em geral, um 
conjunto de regras lógicas que governam a interac- 
ção entre os componentes do sistema. Em certos ca- 
sos, o modelo pode ser explicitado numa relação da 
forma 


onde Y é a variável de saída ou função objectivo 
O Hj Reais - X são as variáveis de entrada ou de 
controle, 

A simulação consiste em introduzir um conjunto 
de combinatórias das variáveis de entrada e procurar 
várias soluções para a função objectivo de acordo 
com as regras lógicas que constituem o modelo. Uti- 
lizando computadores, pode tratar-se um grande nú- 
mero de dados rapidamente e comparar uma série 
de soluções possíveis com a realidade (o que pcde 
levar a um «feed back» sobre o modelo) e intreduzir 
modificações para resolver determinados problemas. 

Segundo o tipo de variáveis de que o sistema 
depende (ou do grau de conhecimento sobre elas), 
a simulação pode ser determinista ou estocástica. 

A simulação estocástica, em que as variáveis de 
que o sistema depende são aleatórias, tem-se reve- 
lado um método útil para resolver uma gama de 
problemas que vão da física nuclear (difusão de par- 
tículas) até à previsão económica. 

Repetindo um certo número de vezes as opera- 
ções lógico-matemáticas do modelo sobre as variá- 
veis de entrada (que obedecem a determinadas leis 
de distribuição de probabilidade), é possível obter 
uma distribuição das resultantes possíveis do sistema 
e sobre esta retirar conclusões que podem ir desde 
uma simples descrição do sistema até à previsão do 
seu comportamento em determinadas condições, opti- 
mizações empíricas, etc. 

Este último caso, em que decisões são tomadas 
a partir da simulação estocástica, é o mais interessan- 
t> e situa-se num campo de inter-fase entre a geoma- 
temática e a investigação operacional. A maioria dos 
exemplos descritos neste trabalho situam-se nesse 
campo. 

Num sistema complexo como é uma unidade mi- 
neira, as primeiras «variáveis de entrada» que temos 
de considerar são VARIAVEIS REGIONALIZADAS 
(dados sobre teores, possanças, etc. em amostras) 
que obedecem a uma lei de distribuição de probabi- 
lidade e a uma função de autocorrelação (covariân- 
cia). 

E a partir deste input básico que um conjunto 
de decisões técnico-econômicas do planeamento mi- 
neiro têm de ser tomadas. 

Se pensarmos em todo o conjunto de operações 
que vai desde a geologia, reconhecimento do jazigo, 
exploração, beneficiação e metalurgia verificamos 
que, ao carácter estocástico particular (*) das Va- 
riáveis Regionalizadas, se vai sobrepor todo um con- 


(*) Tratado em especial na Simulação Condicional. 


junto de interacções com variáveis técnicas e econo- 
micas, algumas das quais são tratáveis por métodos 
determiniísticos, outras por métodos probabilísticos. 

A simulação permite descrever essas interacções 
através de modelos, prever variações no sistema por 
modificações de condições operacionais, comparar al- 
ternativas de funcionamento e optimizar certos pa- 
râmetros. 

É aparente que se trata de uma técnica parti- 
cularmente útil na indústria mineira já que, os cons- 
trangimentos e dificuldades a que esta se tem de 
sujeitar (ao contrário das indústrias transformadoras 
«normais», implantadas num local criteriosamente 
escolhido e tratando, sempre nas mesmas condições, 
matéria-prima de qualidade rigorosamente controla- 
da) tornam difícil o planeamento por métodos em 
que o carácter estocástico não intervenha, 

De facto, a «matéria-prima» da indústria mineira 
encontra-se «in situ», estã sujeita a uma variabili- 
dade importante e nunca é «totalmente» conhecida 
no momento do planeamento. 

Neste trabalho, daremos alguns exemplos do tra- 
tamento deste tipo de problemas através da técnica 
do simulação. 


2. SIMULAÇÃO ESTOCASTICA 
24. Principio de Método 


Seja então Y = X(X, X,...... x ) o modelo que 
estabelece a relação entre a variável de saída Y 
(também chamada função objectivo em certas apli- 
cações da investigação operacional) e as variáveis 
É Ma suubos x. Algumas dessas podem ser Variáveis 
Aleatórias(**), que, em princípio, vamos admitir que 
obedecem a determinadas leis de distribuição de pro- 
babilidade f, (x,), £, (x) ..... f (x). 

O primeiro passo para realizar a simulação es- 
tocástica é fazer uma «amostragem» fictícia das va- 
riáveis X,, X, e... x, isto é, obter 1 realização dessas 
variáveis. 

Se as leis de distribuição de probabilidade dessas 
variáveis fossem uniformes entre O e 1 e visto que, 
nesse caso, a lei de repartição da variável F(x) era 
F(x) =x (ver Fig. 1), números aleatórios unifor- 
memente gerados entre O e 1 seriam os valores da 
realização da variável pretendida. 


ft) Lei de distribuição 


CT. 


Fig. 1 


Distribuição uniforme 


(++) Neste tratamento, vamos admitir que estas variáveis são independentes e obedecem apenas ao constrangimento da 
lei de distribuição de probabilidade. Na simulação condicional tratar-se-á a amostragem de Variáveis Regionalizadas. 
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Verificamos portanto que na base de qualquer 


simulação estocástica está a geração de números 
aleatórios (ver TABELA 1) e 2.2. 

Mas, a partir de um conjunto de números alea- 
tórios uniformemente distribuídos entre 0 e 1 é possi- 
vel obter uma «amostragem» fictícia das variáveis 
XX, cerco obedecendo a leis de distribuição quais- 
quer É (x), É (x,), cmartÊia (34), f(x) através do 
processo seguinte que se descreve na Fig. 2: 


1) Gerar um númro aleatório uniformemente en- 
tre 0 e 1-s BR. 


2) Fazer KR. = F(x) 


3) Obter x=F!(R) 


u 


Fig. 8 


Ficamos assim com n realizações de x,, X, ...... x 
obedecendo às respectivas leis de distribuição de pro- 
babilidade, Introduzindo agora essas realizações no 
modelo, obtém-se uma realização de Y. 

Repetindo esta operação k vezes (ver 2.4 para 
estimar o número de iterações necessárias para cer- 
tos intervalos de confiança) obtém-se uma. distribui- 
ção de Y concretizada num histograma «experimen- 
tal». 

A partir desta distribuição, retiram-se as conclu- 
sões possíveis sobre o comportamento do sistema. 
Estas podem assentar simplemente na Esperança e 
Variância de Y ou noutro estatístico de interesse para 
certos problemas específicos (pode determinar-se a 
probabilidade de Y ser superior a certo limite, por 
exemplo). 

Para esclarecimento mais completo do princípio 
aqui esquematizado, ver [2], [3], [41]. 


2.2. Números aleatórios 


Para obter números aleatórios uniformemente 
distribuídos no intervalo [0,1], podem utilizar-se mé- 


todos mecânicos baseados em qualquer sistema de 
tiragem à sorte, métodos eléctrico-analógicos em que 
um ruído electrónico é convertido num sinal digital 
ou métodos digitais em computador, 

Neste último caso (que é o único prático), os 
números gerados não são verdadeiramente aleatórios 
(designam-se por pseudo-aleatórios) visto que têm 
de ser gerados a partir de um algoritmo qualquer 
(Tabela T). O algoritmo mais utilizado é o das con- 
gruências, em que 


28 E. + c (módulo m) (1) 
isto é, cada número aleatório x, obtém-se do ante- 
rior x, , como o resto da divisão de (a x, , + c) por 
m, sendo a, c e m constantes escolhidas conveniente- 
mente para que o periodo a partir do qual a sequência 
se repete seja suficientemente elevado, 

Por exemplo, tomando 2 como ponto de partida, 
a=3, m=16 e c= 1, a sequência que se obtém 
tem periodo 8, o que é manifestamente insuficiente: 


TR E Z 

x, =— 7 

X, = 5 

X, = 3 

período 8 x, = 10 
x. = 15 

Xe = 14 

x, = 11 

| = 3 


Em alguns computadores, m = ad o Dm ae 
Como os valores obtidos são sempre inferiores a m, 
os números aleatórios pretendidos entre [0,1] obtêm- 
-se simplesmente dividindo por m o output de (1). 

Existem subrotinas para gerar estes números e 
portanto não nos vamos deter neste assunto. Para 
informação competa e para os testes estatísticos que 
garantem a aleatoriedade e independência das sequên- 
cias geradas, ver [2], [8], [4], [5], [6]. 


2.3. Simulação segundo uma lei de distribuição qualquer 


Uma vez obtidos os números pseudo-aleatórios, 
é necessário convertê-los de modo a respeitarem a lei 
de distribuição da variável em causa (*). 

Se o histograma experimental da variável não 
obedecer a nenhuma lei de expressão analítica sim- 
ples, é necessário introduzir, na memória do compu- 
tador, a tabela da distribuição empírica comulativa 
e proceder de acordo com a metodologia descrita 
em 2.1. 

Sempre que possível, esta solução é de evitar 
por razões práticas e económicas e deve interpolar-se 
o histograma experimental através de uma lei de 
distribuição conhecida. 


(*) Este método de «amostragem empirica» é designado por método de MONTE CARLO, 
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Alguns algoriimos para gerar números aleató- 
rios R obedecendo a algumas distribuições mais co- 
muns a partir de Ri; (aleatórios uniformemente dis- 
tribuídos em [0,1]) estão listados seguidamente: 


1) Distribuição uniforme no intervalo [a, b| 
R= Ri(b-a) + a 


-AR 
2) Distribuição exponencial f(R)=)e i ou 
AR 
F(R)=1-€*? 


1 
R = a In (1—-R,) 


3) Distribuição normal reduzida 


Uma aproximação obtida a partir do teorema 
do limite central(*) dá um algoritmo simples para 
obter R 


EK 


Ru, = K /s 


= 
VER 


K é o número de iterações necessárias. 

13 

>» R-—6 o que 
i=1 

dá aproximação satisfatória em muitos casos (lei 
normal reduzida truncada entre —6 e +6). 


Para K = 12, obtém-se R = 


A tabela TI é uma listagem de números aleatô- 
rios segundo a lei normal reduzida. 


4) Lei Normal de média uy e desvio padrão q 


: tz 
R=9 | >, Ry-8 ) +u 


I=1 
EA 
5) Lei de Poisson f(R) ad G=RAO a 


R = K, onde K é o menor inteiro tal que 
K 


1 
2 qn (I-Ry)>1 
1=1 


6) Lei binomial f(R) = (P) p&(1-p)"R 


R=0, 12 .: 
4 O se Ri, <P 
R= > K, 
|=1 onde K = 
E Rs n l se R,>P 


T) Lei de Weibull 


Uma distribuição que se adapta particularmente 
ao tipo de problemas em que se pretende interpolar 
um histograma experimental com uma forma qual- 
quer é a chamada distribuição de Wieibull, cuja lei 
de densidade de probalidade é: 


8 (ty É t=38 0. 
nada co 
e 


E 


e depende dos parâmetros a, 8, y, O que permite con- 
ferir-lhe uma grande variedade de configurações. 

O parâmetro f) (de forma) permite jogar com o 
tipo de distribuição ajustável ao histograma experi- 
mental: 

Para 6=1, a distribuição é exponencial, para 
6 = 3.6 é aproximadamente normal, para 8 = 1.5 é 
aproximadamente lognormal., 

Por escolha conveniente de f, podem obter-se 
distribuições «intermédias» em relação às clássicas 
(Fig. 3). 

O parâmetro wq (de escala) permite comprimir 
ou distender a distribuição com uma certa forma ao 
longo do eixo dos t (a variância aumenta com q). 
longo do eixo dos t (a variância aumenta com q), 
como se vê na Fig. 4. 


til 


O parâmetro y controla a «posição» da distribui- 
ção ao longo da obeissa. 

Um método simples e facilmente programável 
de estimar a, 3 e y descreve-se seguidamente: 


A função de repartição escreve-se 


j id 


Ft) tem a / t>y 


tomando o duplo logaritmo, fica 


a ] =Bin(t-y)-Bim a 


In im | 4 
1—-F'(t) 


(*) Existem outros métodos mais «exactos» de gerar a lei normal, por exemplo, em [3] p. 33 deduz-se um algoritmo 


usando funções trignométricas. 
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Influência do parâmetro & 


fito 
2.5 


2.0+ 


n=0s 


Ls 


Influência do parâmetro ,, 


E + [t-) 
20 3.0 4.0 
Fig. 3 
fazendo Um algoritmo simples de cálculo dos parâmetros 
1 baseado no método dos mínimos quadrados esquema- 
Y= mm) | tiza-se na Fig. 4. 
1—-F'(t) Fig 
X = In (t—-y) ; k 
A=-gha à é um pequeno incremento a dar a y que depende 
Sea da escala dos t e que é ajustável a cada caso con- 
T=EX+A creto. Obtido o y que conduz a um p máximo, pode-se 


INPUT EXPERIMENTAL 


t, Ft) 
CALCULAR 


f, Xn 


En=n-m t+. 


CALCULAR 
an,&n, Zn 
PELO METODO 
DOS MINIMOS 
QUADRADOS 


PRINT *,8,8,2 
Fig. 4 
Diagrama de blocos para ajustar a lei de Weibul 


rodar o programa ainda para obter mais alguns va- 
lores de y na parte descendente da curva p=f(y) e 
obter assim o máximo por interpolação. 


24. Intervalos de confiança 

Introduzidos os valores X, X, ...X simulados 
segundo f,(x,), £,(x,) ..., f (x) no modelo Y=X(x,, 
X,...X ) e fazendo K iterações, obtém-se K realiza- 
ções de Y e o seu histograma experimental, 

É sobre este histograma experimental que as 
decisões vão ser tomadas. 

De uma maneira geral, e na maioria dos casos 
práticos, a distribuição dos Y é normal (de acordo 
com o teorema do limite central). 

Então, os intervalos de confiança a estabelecer 
em Y podem ser construídos a partir das regras da 
estimação dos parâmetros da lei normal a partir de 
uma amostra de dimensão K, 

Por exemplo, suponhamos que as decisões a to- 
mar se referem à esperança da variável de saída Y: 


k 
ça 
Em =T=— bh f 
i=1 
k 
“O der 
E 


(*) Para mais informação sobre o modo de gerar números aleatórios segundo determinadas leis de distribuição, ver 


[71, 181. 


TÉCNICA 451/452 


213 


Admitindo um número elevado de iterações 
(K > 100), os intervalos de confiança para Y cal- 
culam-se por 


onde t, se tira das tabelas da lei normal reduzida 
para determinados níveis de probabilidade 8 (para 
B = 68 %—t, = 1, para 8- 95%, = 1.96, para 
B=99%->t,= 2576). 


Para termos em Y um intervalo de confiança 
fixado, por exemplo E, = 0.05 Y, calcula-se o número 
de iterações necessárias por 


q fp 
f 


Mas este valor mão pode ser imposto à priori 
porque depende de Sº. O processo de cálculo em com- 
putador é, depois de cada iteração, calcular Sº e E 8 
e dar uma instrução de paragem quando E, é menor 
do que um limite pré-fixado que depende do grau 
de confiança necessário para cada problema. E evi- 
dente que se pode estabelecer um balanço entre o 
custo de computação (mais iterações) e o rigor do 
resultado (c, menor). 


3. SUMULA DE APLICAÇÕES 
3.1. Geologia 


Muitos processos geológicos podem ser encara- 
dos do ponto de vista probabilístico, visto que a com- 
plexidade e interdependência entre variáveis dificul- 
tam a representação funcional. 

Por exemplo, as sequências sedimentares podem 
ser modeladas por processos estocásticos, desde que 
s» conheçam as probabilidades de transição de uma 
unidade litógica para a seguinte e as leis de distri- 
buição de probabilidade das espessuras das diferentes 
unidades presentes, 

Se a probabilidade de deposição de uma dada 
unidade S, depende do facto de ter sido precedida 
pela unidade S,, o processo que descreve tal sistema 
é dito de Markov(*). Vamos apenas considerar neste 
exemplo um processo de Markov de 1.º ordem, em 
que o estado do sistema num determinado ponto só 
depende do estado anterior (**). 

Carr et al (1966) [2] estimaram as probabilida- 
des de transição das camadas das Mississipian Ches- 
teran Series por contagem de um grande número de 
casos de litologias que se seguem umas às outras 
numa sequência vertical e construiram a seguinte 
matriz de transição: 


S, S, S, 
arenito Ss 0.1 os 0a 
xisto S, 0.4 0.2 0.4 
calcário S, | 01 8 0.1 


Esta matriz é usada do seguinte modo: se o 
sistema está no estado S, (arenito), a probabilidade 
de surgir de novo outra unidade de arenito é 10%, 
de xisto é 80% e de calcário 10 G.. 

A simulação com base neste modelo de relações 
entre unidades geológicas e leis de distribuição de 
probabilidade pode agora executar-se para prever 
sequências em locais não amostrados, estender as 
observações feitas, desde que haja a garantia de que 
certas propriedades MARKOVIANAS(*) sejam res- 
peitadas. 

O processo de simulação é sintetizado seguida- 
mente: 


1. Escolhe-se um estado inicial para o sistema, 
de acordo com as proporções relativas totais das 3 
unidades litológicas, segundo o esquema seguinte: 


s 
0.85+ * 

Sa 
0.25 

Er 


Se o número aleatório extraído for <0.25 selee- 
ciona-se o arenito com o estado inicial do sistema. 

2. De acordo com a lei de distribuição de pro- 
babilidade das espessuras das camadas de arenito, 
escolhe-se uma espessura para essa camada. 

3. Na 1.º linha da matriz de transição, selec- 
ciona-se a litologia seguinte: 

Se o número aleatório estiver compreendido en- 
tre 01 e 0.9, será um xisto. Da lei de distribuição 
de probabilidade das espessuras das camadas de xis- 
tos, selecciona-se uma espessura para a camada, 

4. Passa-se à 2.º linha da matriz e repete-se 
o processo. 

Na fig. 5 está representado uma sequência pos- 
sivel para a distribuição das camadas. Novas itera- 
ções poderão ser ensaiadas e obtém-se uma «distri- 
buição» de sequências da qual certas decisões pode- 
rão ser retiradas conforme o tipo de problema a 
resolver. 


3.2. Prospecção 


Na prospecção geoquimica, procura-se em geral 
interpretar a distribuição da concentração dos ele- 
mentos pesquisados em águas ou sedimentos e rela- 
cionar as «anormalias» verificadas com mineraliza- 
ções. 


(+) Sobre processos de Markov e sua aplicação na Geologia, ver [17], [18], [19], [20], [21]. 
(++) Stratigraphic Sections, bedding sequences and randon processes — Science, vol. 154 n.º 3753, p. 1162, citado em [211]. 
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Fig. 5 


Simulação de uma sequência sedimentar 
(segundo 21) 


Na maioria dos cases, considera-se «anomalia» 
um «valor elevado» do elemento que se pretende pros- 
pectar. 

Sem discutir por agora esta noção(*), é no en- 
tanto evidente que há um conjunto de factores que 
provocam «falsas» anomalias, cuja investigação no 
campo pode ser evitado se se fizer um tratamento 
preliminar dos dados. 

Tais factores podem ser de vária ordem e são 
em geral considerados subjectivamente pelo geólogo. 

Se for possível quantificá-los e interpretar a con- 
centração dum dado elemento através de um modelo 
do tipo: 

T=f(L HG, CV, M, é) 


(*) 
a teoria das VY. R. 
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onde T é a concentração do elemento 


L é a influência de factores litológicos 

H é a influência de factores hidrológicos 

G é a influência de factores geológicos 

C mede a interferência humana 

V mede a interferência da cobertura vegetal 

M representa os efeitos de corpos mineraliza- 
dos 

e é um termo de erro (ou de factores não 
considerados ) 


Uma técnica de regressão múltipla permite ex- 


plicitar a forma da relação entre T e os diferentes 
factores (que podem não ser estes aqui citados, mas 


dependem do problema concreto). 


O modo de quantificar os factores pode ser o 


seguinte: 


Divide-se a região em estudo num certo número 


de células, como se exemplifica na Fig. 6 (segun- 
do [221]. 


BS. 


9 - 


to 
16 


] É 3 4 5 6 7 8 9 
APS GC 
' 4 ds |, Re | h 
ATI 


1 SA à 


£ a 
papais 


Em cada célula, para quantificar o efeito lito- 
lógico, por exemplo, mede-se a % da célula coberta 
por determinada formação litológica, que passa a 
constituir uma variável do modelo, 

O efeito hidrológico corresponde ao comprimento 
total das linhas de água por célula e o seu número, 
por exemplo. 

O efeito geológico pode medir-se pelo compri- 
mento total de falhas, contactos ou outros acidentes 
estruturais por célula, etc... 

O modelo é então estabelecido pela técnica de 
análise multivariante e um termo de erro aleatório 


Só um conjunto de valores dotados de uma certa «estrutura anómala» deve ser tomado em consideração segundo 
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(calculado através da variância não explicável pelos 
factores controlados e por repetição de amostras) é 
adicionado, 

Simulando, de acordo com o modelo e com a lei 
de distribuição de probabilidade do termo e de erro, 
os valores da concentração nos pontos amostrados, 
ficar-se-á com uma imagem da distribuição que seria 
d: esperar obter nesses pontos se apenas os factores 
considerados tivessem tido influência, 

A diferença entre os valores reais e os valores 
simulados constitui uma medida da qual os facto- 
res perturbantes conhecidos foram retirados e cuja 
distribuição pode fazer ressaltar a existência de pos- 
síveis «anomalias» relacionadas com mineralizações. 

Aplicando o mesmo método(*) numa zona de 
controle para o número ou «valor» de ocorrências 
economicamente interessantes conhecidas, pode simu- 
lar-se uma variável numa zona de estudo (mal 
conhecida) e inferir do seu interesse como objecto 
de prospecção. É este o princípio do <«economical 
appraisal» de grandes regiões. 

Pretendeu-se apenas levantar os problemas que 
podem ser resolvidos por utilização da simulação na 
prospecção. Para informação suplementar, ver [22], 
[23], [24], [25], [26], [27], [28]. 


3.3. Exploração 


3.3.1. Generalidades 


Qualquer método de exploração pode ser dividido 
em operações elementares (subsistemas) com deter- 
minada sequência lógica e interacções previsíveis, 

Assim é possível criar um modelo do sistema 
em que são tomadas em conta as interacções diná- 
micas entre a geometria, geologia e reservas do ja- 
zigo e a alocação do equipamento, métodos e sequên- 
cias de desmonte, 

A técnica de simulação aplicada sobre esse mo- 
delo permite introduzir um método de planeamento, 
estratégias de exploração, e prever modificações no 
output do siztema por variação das condições opera- 
cionais, 

Em qualquer método de exploração, a simulação 
é especialmente aplicável ao subsistema de trans- 
porte dada a necessidade de estudar interferências 
desse processo cíclico com todas as operações ele- 
mentares. 

De facto, em jazigos heterogéneos (zonas «ri- 
cas» e «pobres») a interacção, por exemplo, do teor 
de corte com os custos de transporte (ligados à loca- 
lização do ponto de descarga — rampas, correias 
transportadoras, poços) tem de ser considerada e po- 
dem simular-se diferentes combinações de blocos para 
diferentes estratégias de exploração e várias posi- 
ções dos pontos de descarga, ensaiando uma optimi- 
zação empírica sob certos constrangimentos, 

Por outro lado, se se pretender reproduzir as 
interdependências entre o espectro granulométrico do 
minério desmontado em certas condições (diâmetro 


e disposição dos furos, explosivos, etc.) e o sistema 
de transporte, pode simular-se a influência dos blocos 
de grandes dimensões em estrangulamentos no pro- 
cesso de transporte e prever a dimensão óptima que 
minimize o custo global. 

Assim, justifica-se que se trate o sub-sistema de 
transporte em especial, visto que nele se reflecte, 
como vimos, grande parte das interacções com ou- 
tras operações unitárias. Por outro lado, o problema 
da afectação do equipamento às frentes é um pro- 
blema típico de transporte e difícil de tratar por 
métodos «directos» (sem recorrer à simulação). 

Para informação detalhada da simulação de meé- 
todos de exploração, ver [29], [30], [31], [32], [38], 
[84]. 


3.3.2. Transporte 


É uma aplicação clássica da simulação o estudo 
do sistema de transporte, Uma vez estabelecido o 
grafo de rolagem (em que os pontos do grafo são 
constituídos pelos locais de carga e descarga, bifur- 
cações, etc. e os arcos são constituídos pelas vias) 
é necessário afectar um certo número de transporta- 
dores (camiões, vagonas, correias, etc.) aos carrega- 
dores (scabers, «shovels», pás-carregadoras, etc.) e 
definir uma «disciplina» de tráfego (conjunto de re- 
gras que governam o sistema). 

Na resolução deste tipo de problemas, a simula- 
cão torna-se o método mais adequado, dada a inter- 
acção dinâmica entre os elementos do sistema que 
pode levar à formação de filas de espera (por vezes 
móveis), ao aparecimento de tempos mortos, engar- 
rafamentos, etc. 

A partir da simulação podem formular-se mo- 
deloz que respondam aos objectivos seguintes: 

a) Qual a frota (tipo e quantidade) de trans- 
portadores que conduz a uma produtividade máxima 
(mínimo custo/tonelada), 

b) Qual a organização do sistema («expedição» 
do transportador para pontos de carregamento variá- 
veis ou afectação do transportador sempre ao mesmo 
carregador) que conduz a produtividade máxima. 

c) Qual o perfil de rolagem que, para determi- 
nada frota, é mais aconselhável, 

d) Qual o tipo de equipamento a adquirir para 
determinadas condições, 

e) Qual o efeito da manutenção das vias e do 
equipamento ma produtividade do sistema. 

Suponhamos um grafo de rolagem simplificado 
como se esquematiza na Fig. seguinte 


1 ê 


y - Parque de minério 


ias Escombreira 


Consideremos em primeiro lugar o caso em que 
um certo número de transportadores estã afectado 
sempre ao mesmo ponto de carregamento executando 


(*) Neste caso o que se pretende não é «limpar» a variável de saída de factores interferentes, mas sim «ligá-la» a 


factores que controlem de algum modo o seu surgimento. 
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sempre o mesmo ciclo que se divide nos troços se- 
guintes: 


Tempo de carga; 

Viagem até ao ponto de descarga; 
Descarga; 

Viagem até ao ponto de carga. 


e os Nm 


O tempo comum ao carregador e transportador 
é o tempo de carga e, enquanto o transportador exe- 
cuta o seu ciclo, o carregador retira o minério des- 
montado e prepara-se para nova carga, 

Nestes ciclos existem variáveis aleatórias — tem- 
po de carga e de descarga. O tempo de carga obedece, 
em geral, a uma lei normal ou lognormal, como ze 
exemplifica na Fig. seguinte 


Frequência Relativa 


Tempo minutos) 
Fig. 7 


Tempo de carga (segundo [35]) 


O tempo de descarga é em geral exponencial: 


Frequência Relativa 


Tempo minutos) 
Fig. 8 


Estes histogramas experimentais são obtidos da 
prática do funcionamento do sistema e são válidos 
apenas para aquela combinação transportador-carre- 
gador e para características do minério semelhantes, 
Deve tentar-se ajustar uma lei teórica para poupar 
espaço de memória. 

Quanto aos tempos de viagem, alguns autores (*) 
usam uma simulação determinista baseada na curva 
caracteristica do camião (ou da locomotiva). Esta 
curva permite determinar a força propulsora dispo- 
nivel às rodas em função da velocidade (e de outros 
parâmetros como o peso, declive e resistência à ro- 


(*) Ver por exemplo, [35] e [37]. 
(++) Ver [97]. 
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lagem) e portanto, utilizando um método iterativo 
— num periodo de tempo curto, considera-se a ace- 
leração constante e calcula-se a velocidade terminal 
que dá a força propulsora para o período seguinte 
e consequentemente a aceleração — é possível calcular 
os tempos de viagem para determinadas condições 
de rolagem. 

No entanto, este método conduz sempre a esti- 
mações optimistas do tempo de viagem e é usual 
fazer uma correcção empírica (atribuindo uma certa 
% a mais pelos factores aleatórios) ou colher infor- 
mação na prática que permita estabelecer também 
uma Lei de distribuição de probabilidade para esta 
variável. 

Quanto à tonelagem carregada, esta varia em 
torno da capacidade teórica segundo uma. lei normal. 

O output da simulação será, por exemplo, os 
tempos mortos acumulados para cada unidade e a 
tonelagem total transportada durante um turno, o 
que permite, variando o número de transportadores 
a afectar a cada ponto de carregamento, dimensionar 
a frota para um custo/tonelada minimo e prever 
possiveis modificações na organização do sistema de 
transporte. Quanto a este aspecto, o sistema de «ex- 
pedição» é normalmente mais rentável(**). Este sis- 
tema implica a existência de um ponto de controle 
(e de um sistema de transmissões) que destina o 
transportador ao ponto de carregamento disponível 
ou âquele que o vai estar proximamente. É eviden- 
temente um sistema flexível que se adapta melhor 
à variabilidade dos tempos de carga, à possibilidade 
de avarias nos carregadores e transportadores, etc. 
Há vários tipos de «expedição: que podem ser com- 
parados entre si e com o método «fixo» pela simu- 
lação do sistema, 

No programa de simulação têm de introduzir-se 
algumas restrições ou constrangimentos como sejam 
a velocidade máxima admissível por razões de segu- 
rançca (especialmente em curvas e junções) e regras 
de ultrapassagem e cruzamento, 

Seja qual for o sistema a simular e o output 
pretendido, é sempre necessário introduzir no pro- 
grama um conjunto de regras lógicas que liguem 
entre si as varidveis de estado que são o tempo, as 
coordenadas e as produções. O tempo de cada unidade 
é gerado num «relógio» digital por acréscimos suces- 
sivos segundo as leis de distribuição de probabilidade 
correspondentes a cada operação e esses tempos são 
comparados entre si e com um gerador independente 
(master clock) que fixa o tempo do turno. 

Quanto ao sistema de registo dos tempos mortos 
e produções (contadores) este pode ser a intervalos 
constantes (time slicing) ou só quando se dão acon- 
tecimentos de interesse para o output (evemt sequen- 
cing'). 

Falta agora referir o problema de quantas in- 
teracções são necessárias para uma simulação (quan- 
tos turnos é preciso simular para ter uma certa con- 
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fiança no resultado). De facto, não é possível esta- 
belecer, em geral, intervalos de confiança para simu- 
lações deste tipo (a variedade de situações é de tal 
ordem que não se pode aplicar o método esboçado 
em 2.4.) Neste momento entra de novo a ligação 
com a prática. Não se pode simular um sistema 
completamente desconhecido e portanto, durante a 
colheita de dados, faz-se também a estatística dos 
outputs que se pretendem simular. Antes de fazer 
qualquer modificação no sistema. é necessário «correr» 
o programa nas condições em que os dados foram 
colhidos e verificar quantas iteracções são necessárias 
para que as distribuições do output «estabilizem> e 
se aproximem (estatisticamente) dos histogramas ex- 
perimentais das mesmas variáveis. Esta comparação 
entre as distribuições simuladas e as reais pode ser 
feita por testes estatísticos clássicos (sobre a média, 
variância ou não paramétricos) para saber, com uma 
certa probabilidade, se as duas «amostras» se podem 
considerar provenientes da mesma população. 

Para mais detalhes sobre simulação do transporte, 
ver [35], [36], [37], [38], [39], [40], [41]. 


No método descrito anteriormente, considerou-se 
sempre que os tempos de carga podiam ser tomados 
como realização de uma variável aleatória indepen- 
dente com determinada lei de distribuição de proba- 
bilidade. No entanto, estatísticas disponíveis em algu- 
mas minas (ver, por exemplo, [37]) permitem con- 
cluir que, a tempos elevados, sucedem-se normalmente 
tempos elevados, o mesmo sucedendo a tempos cur- 
tos. Este facto experimental só pode ser interpretado 
pela regionalização espacial do minério que se reflecte 
na «aptidão ao carregamento». Na verdade, duas ope- 
rações sucessivas de carga de minério proveniente 
de pontos vizinhos devem estar ligados por uma 
auto-correlação temporal que é o reflexo da regiona- 
lização espacial medida pelo variograma. Através 
da simulação condicional, este aspecto poderá ser 
tomado em conta. 

Outro problema que não é considerado na simu- 
lação «clássica» das operações de transporte é que, 
por vezes, a variável de saida não é simplesmente 
a tonelagem descarregada (como foi considerado até 
aqui) mas um «indice de qualidade» do minério — 
teor, % de impurezas, etc. 

Neste caso, há que relacionar a «regionalização 
de saída» dessa variável depois de o minério ter 
sofrido todas as operações de arranque, transporte, 
manipulações várias, etc. com a regionalização «in 
situ» dada pela simulação condicional. As operações 
intermédias podem ser tomadas como operadores 
cujas «funções de ponderação» regularizam a variá- 
vel e transformam-na noutra V. R. cuja variabilidade 
é menor e cuja amplitude é maior. Ao simular um 
processo em que o minério vai sofrer homogeneiza- 
ção, este aspecto tem de ser considerado. 
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3.3.3. Planeamento global da exploração 


Em jazidas de minérios complexos, como car- 
bonatitos [43] ou taconites [44] é necessário planear 
rigorosamente a exploração no sentido de pesquisar 
qual a combinatória de blocos que obedece a deter- 
minados constrangimentos sobre índices de qualidade 
limites minimos e máximos em certos componen- 
tes, etc.). 

Estes constrangimentos surgem como especifica- 
ções de entrada na lavaria sobre teores ou outras 
características e é possível, se o número de blocos 
for comportável, simular várias combinatórias possi- 
veis de modo a obter um output de minério desmon- 
tado que obedeça aos constrangimentos. 

Por exemplo, no caso das taconites pode definir- 
-se um índice global de qualidade do minério à en- 
trada da lavaria como 1 ponto no espaço a n dimen- 
sões (uma por cada variável de que depende a quali- 
dade, eventualmente ponderada — por exemplo, teor 
em magmnetite, índice de Bond, rendimento, etc.) e 
definir uma «distância», nesse espaço, de cada bloco 
a esse ponto. Se a mina puder ser dividida em certas 
zonas onde cada pá carregadora se deve manter 
(diminuindo assim as possibilidades de escolha) é 
possivel simular o output proveniente de várias pás 
trabalhando simultaneamente de modo a que a dis- 
tância do índice de qualidade do minério «final» 
(obtido por composição dos diferentes blocos cujo 
fluxo é simultâneo para um ponto de entrada da 
lavaria) ao standard seja minima. Numa determi- 
nada situação, é possível prever qual o bloco a que 
a pá carregadora se deve dirigir no «próximo» turno 
de modo a minimizar tal distância. A aplicabilidade 
do método depende do tempo de cálculo e só em cer- 
tas hipóteses simplificadoras é vantajoso. 


3.4. TRATAMENTO DE MINÉRIOS 


Os processos de tratamento de minérios tem sido 
objecto de um grande número de trabalhos em que 
se pretende estabelecer um modelo que descreva o 
comportamento de certas variáveis, [45], [46], [4%]. 

Sobre este modelo, que relaciona uma variável 
de saída Y com uma série de variáveis de entrada 
ou de controle, é sempre possível executar simulações 
deterministas, prevendo determinados resultados que 
levem à optimização do processo, [48], [49]. 

Não nos detendo sobre este aspecto, vamos con- 
siderar um exemplo em que se executa uma simu- 
lação estocástica, 

Suponhamos que as flutuações previsíveis de al- 
gumas variáveis de entrada se podem prever e que 
o processo está sujeito a uma restrição do tipo: a 
variável de saída (por exemplo, o rendimento ou o 
teor do concentrado) não pode ser inferior a um 
determinado limite —-ou melhor, fixa-se um risco 
para que essa variável seja inferior ao limite pre- 
tendido. 
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Concretizemos num exemplo tratado a partir de 
um modelo construído em [50] por análise factorial: 


Y =26.06-6.29 x,+3.96 x, x,— 6.68 x, x, + 
+ 509 x, x, X, 


onde Y é a concentração de Pb (em %) num pro- 
cesso de flutuação. 


0.05 


0.0 4 


P(Y<45%) 


0.03 


0.02 


0.01 


0.01 | 0.05 


As variáveis «independentes» são: 


A,—0.1 
X, = ———— vindo A, é a concentração em xantato 
0.05 
(lbs/ton ) 
A—1 | 
x, — onde A, é a concentração em OCa 
(Ibs/ton) 
A,—b 
X, = — onde A, é o pH 
1 | 


Admite-se ugora que À, e A, são variáveis alea- 
tórias de lei normal e que se conhece o desvio padrão 
1 Ma 0.05 E (A, 


Ta, = 0.06 E (44) 


O valor máximo de Y obtém-se fazendo x, = —1 
x,=-1 x,=1 e é portanto possível regular o pro- 
cesso para este «óptimo» empírico no domínio de 
validade do modelo. 

Suponhamos que a restrição imposta é a seguin- 
te: Pretende-se que o teor em Pb tenha 2 % de proba- 
bilidade de ser inferior a 45%. Nestas condições, 
visto que a flutuação de A. foi deixada «livre», o 
objectivo é calcular o desvio padrão tal que a res- 
trição seja respeitada para imaginar um sistema de 
controle que assegure que a variável só seja inferior 
a E (A,;) — 2 TA, com a probabilidade de 2.5 % 
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0.10 dp = 02 0.15 


Este é um problema típico, artificialmente sim- 
plificado, em que a simulação estocástica é aconse- 
lhável dado o tipo de restrição imposta e o modelo 
relativamente complexo (devido às interacções que 
impedem o cálculo directo de variâncias). 

O algoritmo utilizado para resolver o problema 
seguinte: Partindo do desvio padrão O (valor fixo de 


Desvio Padrão 


0.20 


Fig. 9 


A.) foram-se dando incrementos sucessivos constan- 
tes e de pequeno passo e para cada incremento foi 
calculada uma distribuição de Y. Em cada distribui- 
ção, leu-se a probabilidade de Y ser inferior a 45 %. 
O resultado destas iterações sucessivas estã esquema- 
tizado na Fig. 9. 


3.5. Análise Económica de um projecto mineiro 


Num projecto mineiro intervém um conjunto de 
variáveis que, no momento da avaliação, não são 
conhecidas completamente. Desde as reservas e teo- 
res que são V. R. cujos estimadores médios podem 
ser calculados pela Geoestatística, até aos preços de 
venda do metal, custos de produção e rendimentos, 
as variáveis econômicas de base são conhecidas com 
determinados intervalos de confiança — trata-se de 
variáveis aleatórias cujas leis de distribuição de pro- 
babilidade, uma vez relacionadas entre si por um 
modelo, podem servir de base a uma simulação esto- 
cástica, 

O modelo a construir relaciona um indicador eco- 
nómico — taxa de desconto interno, lucro, coeficiente 
de rentabilidade, valor acrescentado, etc. — com as 
variáveis económicas de que depende, 

Amostrando as variáveis de entrada pelo Mé- 
todo de Monte Carlo, pode construir-se uma distri- 
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buição do indicador escolhido e, sobre esta, estabele- 
cer, além do valor mais provável, o risco do projecto. 
Assim, associa-se ao parâmetro de decisão um con- 
junto de probabilidades de ele se tornar inferior a 
determinados limites (o risco reflecte a incerteza 
com que são conhecidas as variáveis de entrada no 
momento de tomar a decisão). 

Este método de análise de projecto permite com- 
parar diferentes alternativas, não só do ponto de 
vista do valor esperado do indicador, mas também 
do risco e efectuar análises de sensibilidade à varia- 
ção da esperança e variância dos inputs. 

Às variáveis de entrada tem de ser atribuída 
uma lei de distribuição de probabilidade. 

No caso das variáveis «econômicas» (custos e 
preços de venda) é usual estabelecer uma tendência 
(trend) de evolução previsível no tempo (visto que 
tem de se prever valores futuros) da esperança e 
limites de confiança (em geral o desvio padrão varia 
proporcionalmente à esperança, reflectindo o facto 
de as flutuações previsíveis aumentarem com o au- 
mento do valor médio da variável). 

Alguns autores, [52] utilizam distribuições trian- 
gulares baseadas em probabilidades subjectivas em 
que é atribuída a cada variável um valor modal, um 
minimo e um máximo previsivel. 

Na construção do modelo há que ter o cuidado 
de assegurar a independência entre os diferentes 
inputs (ou estimar as correlações e introduzi-las no 
modelo ). 

Um método em que os valores do teor de corte 
e capacidade anual de produção que conduzem a um 
óptimo da taxa de desconto interno (medida pela 
esperança e diferentes níveis de risco) é descrito 
em [51]. 

Um exemplo de avaliação económico global de 
um projecto mineiro está reportado em [59]. Toma- 
ram-se como variáveis aleatórias o teor do tal-qual, o 
custo de produção, o preço de venda e a taxa de 
desconto (com distribuições triangulares) e como in- 


PROBABILIDADE 
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dice de decisão o coeficiente de rentabilidade defini- 
do como 


= 
my 
DE 


j= 
onde Ro, Do; Lo são, respectivamente, os valores 
actualizados das receitas, despesas e investimentos 
nos N anos do empreendimento. 

A curva de repartição de R está representado 
na Fig. 9 e nela se pode estimar a esperança e desvio 
padrão do indicador (bem como a probabilidade de 
este tomar valores inferiores a certos limites, o que 
mede o risco do projecto). 


R == 
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Bo 247 
2z8 709 
20 344 


93 433 
24 201 
40 610 
T6 493 
61 368 


57 186 
8 253 
64 237 
13 990 
78 822 


89 923 
33 340 
40 881 
89 439 
77 082 


59 593 
10 118 
42 592 
12 059 
29 663 


79335 
82 391 
50 024 
24 892 
83 647 


17157 
92 423 
36 806 
90 033 
35 456 


T3 407 
41 994 
2Zo 298 
20 539 
61 427 


os 021 
19 255 
33 440 
57 546 
21 615 


39 202 
00 822 
do 080 
04 436 
12 171 


11 199 
23 403 
18 623 
83 491 
so 273 


30 500 
oe TTO 
68 711 
29 609 
23 478 


40 218 
14 385 
96 286 
94 400 
5d 382 


31 089 
42 050 
22 222 
28 707 
07 207 


42 582 
33 211 
92 927 
25 701 
05 522 


65 337 
23 287 
90 103 
TS 565 
TO 617 


85 394 
97 118 
84 962 
20 826 
05 014 


do 441 
S7 948 
94 624 
00 387 
TT 938 


80 814 
36 040 
88 461 
15 020 
01 848 


Ta M5 
91 665 
do 606 
27 659 
76 833 


29 170 
09 732 
88 579 
2o 624 
88 435 


73 998 
67 851 
TT 817 
11 062 
34 113 


16 544 
53 763 
02 655 
56 418 
14 598 


20 048 
82 341 
D6 413 
25 815 
31 790 


60 527 
59 466 
45 973 
46 670 
8z2 562 


12 472 
29 529 
39 333 
20 106 
42 941 


11 838 
56 338 
85 207 
69 541 
1 ITA 


31 880 
73 043 
61171 
59 579 
91 936 


51 748 
90 324 
23 306 
09 994 
T6 958 
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TABELA II (segundo [7]) 
Números aleatórios distribuidos segundo uma lei normal reduzida 


1.276 1.218- 0.453 — 0.350 — 0.723 0.676 1.099 — 0.314- 0.914- 0.633 — | 
0.318- 0.799 — 1.664 — 1.391 0.382 0.733 0.653 0.219 0.681 — 1.129 
1.377-— 1.257 — 0.495 0.139 — 0.854 — 0.428 1.322 - 0.315 -— 0.732- 1.348 - 
2.334 0.337 — 1.955- 0.636 —- 1.318- 0.433 — 0.545 0.428 0.297 - 0.276 
1.136- 0.642 3.436 1.667- 0.847 1.173- 0.355 — 0.035 0.359 0.930 
0.414 0.011 — 0.666 1.132 — 0.410- 1.077 — 0.734 1.484 0.340 - 0.789 
0,494 — 0.364 1.237 0.044 — 0.111- 0.210- 0.931 0.616 0.377 — 0.433 — 
1.048 0.037 0.759 0.609 2.043 — 0.290 — 0.404 0.543 — 0.486 0.869 
0.347 2.816 0,464 — 0.632 — 1.614- 0,372 0.074 — 0.916— 1.314 0.038 — 
0.637 0.563 0.107 0.131 1.808 — 1.126 0379 0.610 0.364 — 2.626 
2.176 0.393 0.924 — 1.911 1.040 - 1.168 — 0.485 0.076 0.769 — 1.607 
1.185 — 0.944 — 1.604 — 0.185 0.258 — 0.300 — 0.591 — 0.545 — 0.018 0.485 — 
0.972 1.170 2.682 2.813 1.531 — 0.490 — 2.071 1.444 1.092 — 0.478 
1.210 0.294 0.248 - 0.719 1.103 1.090 0.212 1.185 — 0.338 — 1.134 — 
2,647 0.777 0.450 2,247 1.151 1.676 0.384 1.133 1.393 0.814 
0.398 0.318 0.928 - 2.416 0.936- 1.036 0.024 0.560 — 0.203 0.871 — 
0.846 0.699 — 0.368 — 0.344 0.926 — 0.797 — 1.404 — 1.472 -— 0.118-— 1.456 
0.654 0.955 - 2.907 1.688 0.752 0.434 — 0.746 0.149 0.170 -— 0.479 — 
0.522 0.231 0.619-— 0.265 — 0.419 0.558 0.549 — 0.192 0.334 — 1.373 
1.288 — 0.539 — 0.824 — 0.244 1.070-— 0.010 0.482 0.469 — 0.090 — 1.171 
1.372 1.769 1.057- 1.646 0.481 0.600 — 0.592 — 0.610 0.096 — 1.375- 
0.854 0.535 — 1.607 0.428 0.615- 0.331 0.336- 1.152 - 0.533 0.833 - 
C.148- 1,144 0.913 0.684 1.043 0.554 0.051 — 0.944 — 0.440 — 0.212- 
1.14B- 1.056-— 0.635 0.328 — 1.221 - 0.118 2.045 — 1.977 — 1.133- 0.338 
0.348 0.970 0.017 — 1.217 0.974 — 1.291 - 0.399 — 1.209- 0.248 — 0.480 
0.284 0.458 1.307 1.625 — 0.629 — 0.504 — 0.056 — 0.131 = 0.048 1.879 
1.016- 0.360 0.119- 2.331 1.672 1.053 - 0.840 0.246 — 0.237 1.312—- 
1.603 0.952 — 0.566 — 1.600 0.465 1.951 0.110 0.251 0.116 0.957 — 
0.190— 1.479 0.986 — 1.249 1.934 0.070 1.358 — 1.246 — 0.959 — 1.297 - 
0,722- 0.925 0.783 0.402 — 0.619 1.826 1.272 0.945 — 0.494 0.050 
1.696 — 1.879 0.063 0.132 0.682 0.544 0.417 - 0.666 — 0.104 — 0.253 — 
2.643 — 1.333 1.987 0.668 0.360 1.927 1.183 1.211 1.765 0.035 
0.359 — 0.193 1.023 - 0.222 - 0.616 — 0.060 — 1.319- 0.785 0.430 — 0.298 — 
0.248 0.088 — 1.379- 0.295 0.115- 0.621 — 0.618- 0.209 0.979 0.906 
0.099 — 1.376- 1.047 0,872 — 2.200-— 1.384 1.425 0.812 — 0.748 1.093 - 
0.463 — 1.281 —- 2.514- 0.675 1.145 1.083 0.667 — 0.223 - 1.592 — 1.278- 
0.503 1.434 0.290 0.397 0.897 — 0.973 — 0.120 — 1.594 -— 0.996 — 1.244 - 
0.857 — 0.371- 0.216 0.148 2.106-— 1.453-— 0.686 0.075 — 0.243 -— 0.170 - 
0.122— 1.107 1.039 — 0.636 — 0.860 — 0.895 — 1.458 — 0.539 — 0.159 — 0.420 - 
1.632 0.586 0.468 — 0.386 — 0.354 — 0.203 1.234-— 2.381 0.388 — 0.063 — 
2.072 1.445 - 0.680 - 0.224 0.120 — 1.753 0.571 — 1.933 0.126- 0.034 
0.435 — 0.375- 0.985 - 0.585 — 0.203 — 0.556 — 0.024 0.126 1.250 0.615 — 
0.876 1.227 - 2.647 0.745 — 1.797 1.231 - 0.547 0.634 — 0.836 — 0.719 - 
0.833 1.289 0.022 — 0.431 — 0.582 0.766 0.574 -— 1.153-— 0.520 1.018- 
0.891 — 0.332 0.453 — 1,127 — 2.085 0.722 — 1.508- 0.489 0.496 — 0.025 — 
0.644 0.233- 0.153 -— 1.098 0.757 0.039 — 0.460 — 0.393 2,012 1.356 
0.105 0,171 - 0.110- 1.145- 0.878 0.909 — 0.328 — 1.021 1.613- 1.560 
1.192 - 1.770 0.003 — 0.369 0.052 0.647 1.029 1.526 0.237 1.32B- 
0.042 — 0.553 0.770 0.324 0.489 — 0.367 0.378 0.601 1.996 — 0.738 — 
0.498 1.072 1.567 0.302 1.157 0.720 - 1.403 0.698 0.370- 0.551 — 


224 TÉCNICA 451/452 


NUMERO 451/452 


ANO LIV VOLUME XL 


C. D.U. 519.28:622 


The conditional simulation of regionalized variables 


RESUMO 


A teoria da simulação condicional de variáveis 
regionalizadas é introduzida e as propriedades da 
variável simulada são descritas. 

O método das bandas rotativas para a simulação 
não condicional a três dimensões é apresentado, bem, 
como a simulação a uma dimensão obedecendo a uma 
covariância dada e a um histograma experimental. 

Os problemas de anisotropia, efeito de pepita e 
deriva são introduzidos na teoria e estudadas as mo- 
dificações a realizar para que esses aspectos estrutu- 
rais sejam tomados em conta, 

Alguns exemplos são apresentados para ilustrar 
as aplicações da teoria. 


INTRODUCTION 


Suppose that at the planning stage of a mining 
operation the deposit has been sampled on a given 
grid and from this information it is required to pre- 
diet the fluctuations in mill feed grade so as to de- 
termine whether any blending will be necessary. 
Usually, the limited amount of data available from 
the sampling grid will be insufficient to predict these 
fluctuations. 


The choice of mining equipment may depend on 
the daily fluctuations of the stripping ratio and again 
the available data will be insufficient to predict these 
fluctuations. 


What would be the effects of stockpiling, various 
mining procedures or even different estimation 
methods for determining mining block grades? 

To answer these questions would require knowle- 
dge of the measured variable at all points in the 
deposit. As this is obviously impossible, an alterna- 
tive would be to simulate the measured variable at 
all points in the deposit in such a way that it con- 
forms to the known characteristics of the variable 
as deduced from the available data. 

However, it must always be born in mind that 
simulation is not estimation — its purpose is to pro- 
vide unbiased variances, 
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SUMMARY 


The conditional simulation theory of Regionali- 
zed Variables is introduced and the basic properties 
of the simulated variable are described. 

The turning bands method for non-conditional 
simulation in a three-dimensional space is presented 
as well as the one-dimensional simulation generated 
by a given covariance function and experimental his- 
togram. 

Problems of anisotropy, nugget effect and drift 
are introduced in the theory and modifications for 
taking these structural features into account are 
done accordingly. 

Some case stories are presented for application 
of the theory. 


These notes are based on Journel (1974) (refe- 
rence number [1]) and this source quoted freely in 
the text. 


THE PROPERTIES OF A CONDITIONAL SIMULATION 


The available information about the deposit con- 
aists of: 


1 —'The values z(x,) at each of the sampling 
locations. 

2 —'The spatial autocorrelations between data 
Z(x,) and z(x,). These autocorrelations are 
expressed in the form of a covariance or 
variogram function, 

3 — 'The experimental histogram of the data va- 
lues, 


Thus, the simulated variable, z (x) must: 


1 — Be equal to the experimental value at each 
of the data locations, i.e, ZÃx,) = 2(x,), 
V x, € 8. In other words, the simulation is 
conditional to the data; 

2-— have the same spatial autocorrelation as the 
data values i.e. the same covariance or va- 
riogram; 


(*) B. Sc. M. Sc. A, M. I. M. M., M. Can. 1. M. M. Ph. D. 
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3 — have the same histogram as that of the data 
values, 


THE THEORY OF CONDITIONAL SIMULATION 


At each point x in the deposit (D) the neighbour- 
ing data values, z(x,), can be used to form a Kri- 


ged estimate, z 4X), of the true value z(x): 


“ à 
Z, AX) = > A, Z(x,) 


| 


The difference between the kriged value and the 
true value is the kriging error: 


e, (x) = 2(x) — Z 4x) 
Z(x) = 2,4%) + e (x) 


Now, if all these values are considered as reali- 


sations of random functions, Z(x), Z (x), E, (x), 
then the kriging error is orthogonal to the kriged 


estimate, i.e. E Z, (x) E, (x') | =0. 


Now consider realisations, z (x), of another ran- 
dom function, A 4x), which has the same expectation 
and covariance as Z(x) but is independent of it. 
& (x) could result from non-conditional simulations 
of Z(x) i.e. a simulation that simply ensured the 
same expectation and covariance. 

Using the values of z (x) at each of the same 


data locations x, € S, each location x € D can be 


kriged to give an estimated value, z ' (x): 


2.41%) = ke à 2, (x) 


where the weighis à are the same as the previous 
set of weights because the same data pattern is 
being used. 

This leads to the following expression: 


Z0)=2Z,0)+E (x) 


where the kKriging error, E, (x) = Z (X)-Z (x) 


has the same expectation and covariance as E (x) 
and thus the random function: 


Z (x) =2 (x) +E, (x) 


k 


has the same expectation and covariance as Z(x). 
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This last random function is obtained by adding the 


kriging error E, (2) — [2,00-25, 00 | to the 


kriged estimate 7 (x) of the initial random function 
Z(x). These results are due to the orthogonality 
property mentioned above. 

We now have a method for conditional simula- 
tion which can be summarised as follows; 


1 — Carry out a non-conditional simulation of 
the regionalised variable for all locations 
x € D. This results in a set of values (Z (x), 
y x € D). These simulated values have the 
same expectation and covariance as the data 
values [2(x,), Y x, € 8). 


2 — Using the data values z(x,) at the sampling 
locations x, ES Krige all locations x € D. 
This provides a second set of values 


(4 (x), X€ D4. Remember from the exact 
interpolation property of kriging that 


2,4X,) = 2(x,) v x, € Ss 
9 — Using the simulated values z (x,) at the 
sampling locations x, € S, krige all locations 
x€ D. This provides a third set of values 
(2 (x), x € D). Remember from the exact 
B,k 


interpolation property of kriging that 


Z y (x)=2z (x) V x, é Ss 


Then the realisations: 


z (x) A o [2,00-2,, 00] 


are conditional to the real experimental data since 
at any data location x E S: 
21x) = 2. 4X) and 


+ Z (x) = z(x, ) 


and therefore: 


ZAX,) = 2(x,) V x, €S 


In addition, the corresponding random functions 
£ ix) and Z(x) have the same expectation and cova- 
riance, 


NON-CONDITIONAL SIMULATION: THE TURNING 
BANDS METHOD 


This section deals with the generation of a reali- 
sation z (x), of the random function Z(x) with a 
given covariance, 

The turning bands method consists of reducing 
any three-dimensional simulation to simple one-di- 
mensional simulations on lines that turn in the three- 
-Cimensional space. The method was developed by 
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Matheron and colleagues at the Centre de Morpho- 
logie Mathématique de Fontainebleau [1], [2], [9]. 


— Consider the line L, defined by the unit vector 
u, in the three-dimensional space as ehown in the 
following figure: 


On line L, consider a one dimensional random 
function Y(u,) which is second order stationary with 


e, ? yqtua) ytua 
yu =. 
X X 
1 2 3 
ytuo) 
q 
15 
X X X 


/ 


E) 
of vector h on these axes, 


The three-dimensional random function Z(u, v, w) 
defined by 


z(u, V, w) = y(u) 


is also second order stationary with zero expected 
value and a three-dimensional covariance of: 


OC, GBah,h)=0th) 


In practice, to generate a realisation z of Z, the 
value of the one-dimensional realisation v(u) at the 
point u on line L, is assigned to some given small 
Interval centred on the point u. This interval is then 
projected at right angles to L, into the three-dimen- 
slonal space. Every grid point that falls within this 
Interval, or band, then receives the value y(u). This 
is illustrated in two dimensions in, Fig. 1. 

Suppose that we wish to simulate a given regio- 
nalised variable at each of the prid points marked X. 

On line L we generate the values of y,(u,) at 
each of the points (e). These points are regularly 
spaced at an interval b. 'These intervals are then 
projected at right angles into the two dimensional 
space. Grid point 1 receives the value of y (u,), grid 
points 2,7 and 13 receive the value y (u,) and so on, 

In the three-dimensional space the intervals are 
projected by means of planes as shown in. Fig. 2. 


£ 
X X X 
d 5 6 
X X X 
10 n 1 


Fig. 1 


a zero expected value E(Y)=0 and has a one-dimen- 
sional stationary covariance C, (h ). 
l 1 


= > —s> 
Let (uu, v, w ) represent a unit system of 


orthogonal axes and (h,, h, h,) be the projection 
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All grid points that fall within the band centred on 
the plane (u, = constant) receive the value y(u,). 


We then consider N directions u,, u, ... u, uni- 


formly distributed on the unit sphere. A new realisa- 
tion y(u,) of a random function Y,, equivalent to Y, 
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y(ui;) 


Fig. 2 
is generated on every line L,, the N random functions 
(Y, i=1, N) being independent. For every one- 
“dimensional realisation y(u,) there is a correspon- 
ding three-dimensional realisation Z, (U,, V,; W,) = 
= v(u,). 


We then consider the sum 


which is a realisation of a three-dimensional random 
function Z(x, x,, X,) which is second order stationary 
with a zero expected value and a covariance which 
tends towards the following isotropic covariance as 
the number of lines, N, tends towards infinity: 


. —» — 
C,(h, h, h,) = €(h) = | C(<h, u>) du 


1 unitsphere 


—, —s 
where < h, u>=ãh denotes the projection of vec- 


tor h on axis u and n=/ hº += há + hº is the modu- 
lus of vector h. 

The above integral can be written in terms of 
spherical co-ordinates as: 


2% / 
C(h) = | 


T/a 
do | C, (|hcos q|) sin q dy 
0 


0 


which can be inverted to give: 


1 
C, (9) = —— 


s C(s) 
2x Os 
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In practice N is never infinite. The method 
adopted at the Centre de Morphologie Mathématique 
de Fontainebleau is to take the lines that join the 
mid-points of opposite edges of a regular icosahedron. 

The sum 


will then have the covariance: 


15 
> —s 
em= Sec RBT>) 
j=1 
which will differ from the limit covariance 
2x h 
Sa 4 eins À C, (s) ds 


The approximate relationship between these two 
covariances is 


15 
C(h) = —— O(h) 
24 


In summary, to generate realisations Z(X,, X,s X,) 
of a random function Z(x,, x,, x,) with an isotropic 
three-dimensional covariance function C(h): 


(1) — Find the corresponding one-dimensional co- 
variace function C, (h). 


(2) — generate realisations, y, (u,), of the one- 
-dimensional random function Y, (u,) 
which has a covariance C (S) on each of 
the (L,, 1=1, 15) lines. 


(3) — project these realisations into the three- 
“dimensional space and sum them on the 
simulation grid locations. 


The only remaining problem is the determination 
of C (s) and the generation of the realisations of the 
one-dimensional random functions Y (u) where 


C,(s) = E (X(u) Y (u+ 8)) 


E (Y(u)) = O 


ONE DIMENSIONAL SIMULATIONS 


The most general method for generating these 
1-D realisations is by means of spectral analysis but 
in practice this may require very difficult Fourier 
transforms. Simple methods exist for the particular 
cases in which the covariance C (s) can be written 
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as a convolution product of a known function f(u) 
by its transpose 


[2 
Glm= £* f = f f(u) f(u + s) du 
— 0 


Then a moving average method with weighting 
function f(u) can be used to generate the Y(u). We 
start with a stationary 1—-D random measure T(dv). 
Then for each measure T(dv) there is an associated 
1-D random function Y(u) defined by: 


[e 
T(u) = f fiu + v) Tí(dv) 
— DO 


is a stationary random function with covariance 
C,(s) = É * f 
SOME EXAMPLES 


(1) The spherical covariance 


The three-dimensional spherical covariance func- 
tion is written as 


ir dh hº wa 
C(h) = “2a Li 2 u3 “ 
() h> a 
now in the one dimensional space: 
, 1 ) 
G, (s) = Ds Vs s C(s) =" 
1 à õs s8 
2x Õs i ( 2a gs) 
fo SB, 2 s<a 
a 2x a as 
0 5s>a 


This covariance can be expressed as the convo- 
lution product: 


Cutp=2*4 
where the weighting function is: 
V BE a 
f (u) mas 


O otherwise 


<u< a, 


EPs af 
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demonstration: 


GC (s)=[s [= t(u) f(u-+s) du = 


nes n's + & 
dra 
fu (u+s) du= 


5 q 
o VS 
mas “ -Alg+B 
" 3] us a 
<S|d 2 Lata 


» K 1 E J 58 
- q T giz a a as 
6 K 


Thus, neglecting the multiplicative factor Ve 
the 1-D random function, Y(u), for the spherical 
covariance function is expressed as: 


u-+ als 
Y(u)= | (u—y) T (dy) 


In practice, this procedure consists of assigning 
a value, Z,, drawn at random from a uniform distri- 
bution, to each of a set of points x, marked at regu- 
lar intervals along the line. Then to each point x, 
the value 


do Vihsi 


k=— als 


is assigned. 


(2) The exponential covariance 


1 0 
now. C, (5) = 2. Ea s C(s) = 
1 0 -—as 
= 2x dOs pera is 
E (1—as)e * VV 020 


This covariance cam be expressed as the convo- 
lution product: 


Cl) =f"*1 


where the weighting function is: 


yZÉs sas 
f(u) = do 
0 
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In practice, this procedure consists of assigning 
a value Z,, drawn at random from a uniform distri- 
bution, to each of a set of points x marked at regular 
intervals along the line. Then to each point x the 
value 


a 
nd = 
y, = > (ake E as 
k — U 
is assigned. 


GEOMETRIC ANISOTROPY 


Thus far we have dealt with covariance functions 
which are isotropic in the three-dimensional space. 
In fact, the derivation of the one-dimensional cova- 
riance function from the three-dimensional covariance 
function assumes isotropy and the derivation is 
practically impossible unless this assumption is made. 
However, in practice a reasonable approximation can 
usually be made by calculating the ranges (2, 
1=1,15) corresponding to each of the (L,, i=1,15) 
lines. 

An alternative is to simulate an isotropic cova- 
riance function on a given grid and then multiply 
this grid by the ratio of anisotropy. In two dimen- 
sions (for the sake of a simple example) if the range 
is twice as long in the E/W direction than in the 
N/8S direction then the co-ordinate system would be 
transformed by multiplying the E/W grid spacing 
by two. 


THE NUGGET VARIANCE 


The nugget constant appears as a random com- 
ponent of the total variance — it represents a micros- 
tructure of very small range i.e. practically no spatial 
auto-correlation. 

It can be simulated by drawing numbers from a 
uniform distribution with a variance equal to the 
nugget constant, A uniform distribution on the inter- 


wx* 
val [o, x] has a variance equal to -—-—. Therefore, 


1 
if the nugget constant is C we draw numbers from 
a uniform distribution on the interval [o, y 12C]. 


ZONAL ANISOTROPY 


If, for example, the variogram function consists 
of an isotropic component, Y, th), in the three-dimen- 


stonal space and a second component, y (vh2 + h2), 


depending only on the horizontal distance, then these 
two components can be simulated independently and 
the resulting realisations added together. 


SPECIFYING A HISTOGRAM 


Thus far we have satisfied the first two condi- 
tions of the conditional simulation process l.e. condi- 
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tionality, and covariance. The final condition is to 
ensure that the histogram of the simulated values is 
the same as that of the data values. 

As we are summing a large number of values 
on each line the resulting realisations, z(x), will tend 
to follow a normal distribution because of the central 
limit theorem. 

In general, to impose a given histogram on the 
simulated values we must find some transform, q, 
which will transform the data values to values which 
follow a normal distribution: 


v(x) = g[2(X))] Vxes 
where y (x) follow a standard normal distribution. 


We then determine the covariance function cor- 
responding to the y(x,) and simulate these values. 
When the simulation is completed we take the in- 
verse transform: 

ZAx) = q [Y. (X)] yY xe D 

The values z (x) will then have the required his- 
togram and covariance function. 

The transform, &, can be achieved in a number 
of ways. When a known distribution cam be fitted 
to the data, the transform can be simply determined. 
For example, if it is known that the data follow a 
lognormal law then the transform is achieved by ta- 
king logarithms: 

v(x) = log[2(x,) ] V xos 

We determine the covariance function of the 
y(x,) and simulate the values, y (x), xED. When 
the simulation has been completed we take the in- 
verse transform: 

24x) = SAX) Y xeD 

The values z (x) then follow a log-normal distri- 
bution and have the required covariance. 

When an exact distribution function cannot be 
fitted to the experimental histogram, a graphical 
transform using the familiar Monte Carlo techniques 
can be employed. 

The most general method is to use hermite poly- 
nomials to obtain a numerical transform. The her- 
mite polynomials are defined by: 


x? /92 dº Miaidadd 


Hix)= 8 
n dx” 


Now if the random function Z(x) is strictly sta- 
tionary and has a distribution function F(z) then a 
transform, q and its inverse, 4", can always be 
found such that 


Z(x) = é [Y(x)] and Y(x) = 4 [Z(x)] 
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where the random function Y(x) has a standard nor- 
mal distribution, 

The transform q can be expressed in terms of 
hermite polynomials as 


so 


3 


=0 


g(y) = 


H (y) 


where the weights v are given by: 


o 
[ +) H (yig(y) dy 


— (O) 


Fa = 


and g(y) is the standard normal density function 
and y = yíx) 

Now Z(x) = g[Y(x)] implies that Fiz) = G(y) 
where F(z) and G(y) are the cumulative distribution 
functions of Z and Y respectively. The values of the 
transform cam be determined for W x, € S, from: 


z(x,) = gly(x,)] 


where y(x,) are defined by 
F[z(x,)] = Gly(x,)] 


Knowing the pairs of values [z(x,), gly(x,)]]. 
It is a simple matter to use any standard numerical 
integration method to obtain estimates of the weights 
“, The transform q(y) is expanded up to order N 
where N depends on the shape of the distribution 
(experimental histogram) of the z(x,). Terms up to 
and including v, and large sample sizes are needed 
(and will almost always provide) a good approxima- 
tion to g(y). Fewer terms and small samples can lead 
to a poor approximation, including negative proba- 
bilities. 

The hermite expansion of q(y) should be such 
that: 


where m and c' are the mean and varlance of the 
u(x,) V x, €5. 

Once the function s(y) has been determined as 
a hermite expansion, the equation z(x,) = gly(x,)] 
ca be solved mumerically or graphically for each 


z(x,) to determine the inverse transform, & (2). 


SIMULATING VARIOUS VOLUMES 


Usually the available structural information will 
be on a particular support such as drill cores. Using 
this information we may wish to simulate drill cores 
of various lengths, blocks of various sizes and other 


TÉCNICA 451/452 


types of samples, One method of doing this is to 
deduce the structural characteristics for a point 
support and simulate the variable on a point support. 
Various volumes can then be simulated by averaging 
the values on all points which fall within these volu- 
mes 

A more direct way is to assume invariance of 
the type of distribution for various lot sizes. For 
example, if the point support distribution is lognor- 
mal then the block-support distribution will also be 
lognormal with the same mean but a smaller varian- 
ce. This is known as the uid (or invarlant) 
law hypothesis [14]. 

As before it is assumed that the sample grade 
values z(x) can be expressed as standard normal 
transforma: 


z(x) = gly(x)] 


which can be expressed in terms of hermite polyno- 
mials as: 


o 
sm= SL Et) 
dito. 4 


Suppose that the required block grades are z(v), 
vt D, The invariant distribution hypothesis amounts 
to assuming that the standard normal transform, q,, 
of the z(v): 


z(v) = 4, [wlv)] 


has the same hermite polynomial expansion as q 
except for a correction factor: 


So ae 
> Tu tn Ho (wW) 


n=0 ds 


BW) = 


The c, can be determined from the variance 
relation of the hermite polynomials [11]. 

The most general way of estimating block grade 
distribution is by means of transfer functions [10], 
[11], [12]. The (direct) transfer function of blocks 
v, in the deposit (or zone) D, where D = Uv, is the 
conditional distribution of the block grades, z( V)), 
given the sample grades z(x,), x,€ 5. 

Determination of a transfer function consists of 
forming a disjunctive kriging estimator of the pro- 
bability density function of the block grades given 
the sample grades. A full account of disjunctive kri- 
ging (dk.) can be found in [10] and [18]. 

The disjunctive kKriging estimator of the transfer 
function (or the psuedo-density function) of the 
blocks V, ig 


Árti 
fax (5) = 8 (5) PR in > 4 E Ent] 
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where g(s) is the standard normal density function 


and the coefficients Nú are determined by solving the 
following system n = 1, 2,... N times: 


l n BA ecra: 
> A po Ey(xD, VOX = [ytxo), wív)] v k 


where p(h) is the covariance function of the standard 
normal variable y(x) and 


p" [y(x); y(x0] = Lolx, — x97º 


r(h) is the covariance function between the standard 
normal variables w(v) (blocks) and vy(x) (samples) 


and rº [yíx), w(v)] is the mean value of r (X,—V,) 
v v, ED. 


DRIFT 


Sometimes the spatial variation of the simulated 
variable include a directional trend due, for ins- 
tance, to progressive impoverishment or enrichment 
towards the borders of the deposit, In this case, the 
random function Z(x) cam be considered as the sum 
of two random functions: 


Z(x) = T(x) + R(x) 


where T(x) represents the trend and R(x) is the 
residual fluctuation around the trend. 


v(h 


03 


020.206 


o2 | a 


01 
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Realisations of the trend and its residuals are 
simulated independently, see, for instance, [15]. 


EXAMPLE 1 


This first simulation [6] was carried out to 
forecast departures from planning in open pit design 
and grade control. 


Using the known characteristics of a porphyry 
copper deposit a simulation of the copper grades on 
point supports was carried out. 


The variogram of the point grades was isotropic 
in the horizontal plane with a range of 400 ft while 
the vertical range was 100 ft. There was a nugget 
constant of 0.05 (%)” Cu and a total variance of 0.20 
(%)* Cu The mean grade was 1.19 % Cu. 


A total of 67,200 points were simulated on a 
0 ft X 50 ft xXx 12 ft grid over a field measuring 
2000 ft X 2000 ft X 504 ft. The variograms resulting 
from the simulation are shown in, Fig. 3. 

These simulated points were averaged along 50 
ft lengths to simulate drill cores on various drilling 
grids. Points within blocks measuring 100 ft X 100 
ft X 100 ft were averaged to provide one possible 
set of «real» block values. 

Various estimation methods, such as inverse - 
squared distances and kriging, were then used to 
estimate the blocks and calculate grade/tonnage 
curves. These estimates were then compared with 
the <real> values. 


500 600 700 800 500 | 


Fig. 3 
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The blocks were then kriged on different drilling 
grids and the resulting estimates were used to deter- 
mine optimal open pits. 

The full study can be seen in reference [6] 


EXAMPLE 2 


This example is taken from Jourmel [1| and 
concerns the Prony nickel deposit in New Caledonia. 
A total of 176 cores from boreholes on a 100m 
grid were available. This data was to be used to 


%k histogram 
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predict flutuations in grades and stripping ratios. 
As this information was insufficient a conditional 
simulation was carried out on a 10m grid providing 
a total of 17,600 dril cores. 

The variables simulated were 


overburden thickness 

mineable thickness 

Nickel metal accumulation over the mineable 
thickness. 


A summary of the main results is shown in the 
following figure. 


m = 23,08 
Experimental data 
D* 
= =— 0,08 
m 
m =23,7 


— — — — Conditional simulations 


meters 


so 50 


Theoretical model 


----— Conditional simulation 
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From Journel [11]. 
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SUMÁRIO 


Dado o carácter prático deste trabalho, e estando 
ele integrado no Seminário de Geomatemática, onde 
outros trabalhos apresentam os conceitos teóricos da 
simulação, estes aspectos não são incluídos aqui. 

Por outro lado, sendo este processo de simula- 
ção uma aplicação à Análise Económica de Projectos 
admite-se que: 


— São conhecidos os conceitos necessários 
ao Discounted-Cash-Flow. 

— Um projecto /exemplo está disponível e já 
foi explicado. 

— São conhecidos os princípios básicos da 
Informática. 


Para satisfazer os dois primeiros pontos o leitor 
deve remeter-se qo anexo, onde de uma maneira su- 
cinta se relembram os conceitos, se introduz o pro- 
jecto/exemplo, e se analisa o efeito de variações dis- 
cretas (Análise de Sensibilidade). 

Neste trabalho propriamente dito, expõe-se a es- 
trutura do processo de simulação, apresentando o 
projecto /exemplo na óptica dos conceitos probabilis- 
ticos (Análise do Risco). 
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HEITOR SIMÕES SANTOS * 


SUMMARY 


As this paper intends to be a practical one as 
well as a part of the present semminar, where theo- 
retical concepts on simulation were already intro- 
duced, the former part was skipped. 

On the other hand, being «This Simulation Pro- 
cess» concerned with the Economic Analysis of the 
Projects three points are assumed: 


— Modern concepts of Econ. An. of Proj. 
(DCF) are known. 

— A project/example is available and al- 
ready explained. 

— Basic concepts on the computer science 
are required. 


Excluding the last point, the reader is reverted 
to the Annex, where a short introduction on DCF 
is available (as a reminder) as well as the Project/ 
(Example being introduced in order to clear the 
structure and the concepts explored on discrete va- 
riations (Sensitivity Analysis). 

Therefore, the present paper is concerned with 
the structure of a simulation process, presenting The 
Project/Example under the probabilistic concepts 
(Risk Analysis). 


1 — INTRODUÇÃO 


Integrado neste seminário, o presente trabalho 
tem objectivos limitados que o confinam à «prática 
do processo de simulação» aplicado à análise de in- 
vestimentos. 

Trata-se, portanto, de um trabalho que, reque- 
rendo o conhecimento e a capacidade de utilização do 
Cálculo Automático, se apoia fundamentalmente nos 
conhecimentos de Análise de Investimentos e dos 
Processos de Simulação. 


Abordaremos assim, sempre sucintamente e com 
objectivos específicos, os temas propostos, fazendo-os 
seguir da descrição e elementos/base do exemplo, 
para depois apresentar os resultados e as conclusões, 


* Engenheiro de Minas (I. 8. T.), M. Sc. (University of London), Eq. Prof. Aux. (Universidade de Coimbra). 
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2 — PRINCÍPIOS DE ANÁLISE ECONÓMICA DE INVES- 
TIMENTOS (D. C. F,) 


Este tema é pressuposto ser do conhecimento dos 
participantes, mas uma exposição sucinta encontra-se 
nas primeiras quatro páginas do ANEXO I. 

A técnica aqui seguida é portanto a do «Dis- 
counted Cash-Flow» produzindo como indicador prin- 
cipal a Taxa de Desconto Interna (TDI = «rate of 
return») e como indicadores complementares o Má- 
ximo Capital Investido (MCI) e Tempo de Recupe- 
ração do Capital (TRC). 


3 — PROCESSO DE SIMULAÇÃO (MONTE-CARLO) 


Perante a «Introdução à Simulação» já apre- 
sentada, tendo em atenção os Cap. 1, 2, 3 e muito 
especialmente o sub-capítulo 3.5, apenas especifica- 
remos que o processo aqui seguido se baseia nas 
técnicas de  Monte-Carlo, sendo portanto um pro- 
cesso de simulação estocástico. 

Quanto às distribuições adoptadas para as variá- 
veis aleatórias, são do tipo triangular. 

Mas, dado o carácter prático deste trabalho, 
tornaremos mais evidente a estrutura de um pro- 
cesso de simulação. 

Na realidade podemos considerar três fases dis- 
tintas: 


— Geração de números pseudo-aleatórios de acor- 
do com as distribuições definidas para cada 
variável, 

— Modelo matemático que traduz o método de 
análise em questão, produzindo como resulta- 
dos os chamados «indicadores». 

— Análise estatística dos conjuntos de indicado- 
res calculados. 


Este esquema será pormenorizado e evidenciado 
no Cap.º seguinte. 


GERAÇÃO DE 


| NÚMEROS ALEATÓRIOS R MODELO 
HATEMÁTICO 


DE DISTRIBUIÇÃO 


APROPRIADA 


ela não tem outro objectivo que não seja a identi- 
ficação e situação do presente caso. 

Os programas foram escritos em BASIC e estão 
formulados de acordo com as fases do processo de 
simulação, a que já fiz referência, encadeando-se 
contudo automaticamente, 

A Fig. 14 apresenta a estrutura do processo em 
termos informáticos. 

As razões que levaram a particularizar o pe- 
queno esquema do cap.” 3 do modo apresentado não 
estão no âmbito deste trabalho mas tem conse- 
quências de carácter prático das quais ressalta a 
possibilidade de, em caso de acidente operativo, ou 
de se pretenderem alterar os elementos/base do 
modelo, não ser necessário reiniciar todo o processo, 

Estas vantagens não serão de desprezar se, 
mesmo quando se trabalha com computadores, re- 
cordarmos que o modelo terá de ser «corrido» um 
número elevado de vezes — que pode atingir as cen- 
nas — de acordo com a precisão desejada. 


5 — ELEMENTOS/BASE DO PROJECTO/EXEMPLO 


O projecto/exemplo sobre o qual fizemos incidir 
este estudo encontra-se descrito no ANEXO I, com 
início na 5.º página. 

Do tratamento, que ainda no ANEÃO 1, mas 
nas páginas seguintes à descrição, o projecto /exem- 
plo é submetido, só nos interessa concretamente, 
aquele que se encontra da 5.º à 8.º página e que é 
chamado de «cash-flow básico», 

Convém contudo referir, embora isso se encon- 
tre em evidência na 9." página, que daí para diante 
se procedem a estudos de sensibilidade, que, como 
se sabe, são estudos relativos a variações discretas 
dos elementos/base do projecto/exemplo, aqui tra- 
duzido pelo modelo matemático. 

Assinalâmos porém este facto pelo que ele re- 
presenta como conceito diverso do que preside ao 
processo de simulação (distribuição probabilística). 


ANÁLISE 
ESTATÍSTICA 
DOS CONJUNTOS DE 
INDICADORES 


FIG. 1.3 - FASES DE UM PROCESSO DE SIMULAÇÃO 


4 — CÁLCULO AUTOMÁTICO 


A necessidade da utilização de computadores 
aparece como evidente, pelo menos como único meio 
prático de realização dos processos de simulação. 

Por isso, esta pequena nota se justifica, mas 
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5.1 — As variáveis aleatórias 


Das variáveis que enformam o modelo do pro- 


jecto/exemplo, nem todas se podem classificar de 
aleatórias, pelo menos se pretendermos encarar o 
problema em termos práticos, 
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| PROGRAMAS PROGEAMAS PRINCIPAIS FICHEIROS 
AUXILIARES 
E) Dera 
DATA - Limites dos coeficientes 
a das variáveis aleatórias 
E FUNÇÃO - Gera num. aleatórios de 
igual distribuiçao,trans- 
forma-os na dist. trian- 
gqular respectiva e escre- 
ve-os no ficheiro SDCF12 
SD era 
DATA »- Elementos/base do modelo 
cá FUNÇÃO - Afecta os elem./base com 
H os coef. respectivos 
E) - Calcula o modelo definin- 
pe do o cash-flow 
E | a o - Calcula os indicadores e 
Ed escreve-os em ficheiros STRC23 
TRC SMCT23 
| S DCF 3 
=” AR A A = Am 
ni | DATA - Contida nos ficheiros 
A FUNÇÃO - Analisa a distribuiçao 
Ê dos indicadores e dese- 
| nha o gráfico respectivo 
ES To , 
ma DIST. CUM. DIST. CUM. DIST.CUM. | 
ES TDI MCI TRC 
O! 
| 


FIG. 1.4 - ESTRUTURA INFORMÁTICA DO PROCESSO DE SIMULAÇÃO 


Assim admitimos que apenas cinco dessas va- 
riáveis justificavam um tratamento aleatório, ou me- 
lhor, para ser mais preciso, cinco conjuntos de va- 
riáveis. 


Elas são: 1) — Preço de venda dos pehble 
2) — Preço de venda dos smalls 
3) — Investimento em equipamento 
4) — Capital operacional fixo 
5) — Investimento em construção civil 
6) — Custos operacionais 
7) — Capital operacional móvel 


onde 4) se encontra directamente ligado a 3), e 7) 
se encontra directamente ligado a 6), constituindo 
assim os cinco conjuntos. 

Os valores mais prováveis das variáveis são os 
que deram forma ao modelo base. 
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O problema que se põe inicialmente é o da 
definicão do rigor com que esses valores foram cal- 
culados e consequentemente a atribuição a cada uma 
dessas variáveis, de coeficientes correctivos determi- 
nados aleatoriamente e com distribuições apropria- 
das (triangulares), cujos limites inferior e superior 
hã que definir. 

Assim, teremos para os cinco grupos de va- 
riáveis: (Fig. 1.5.1). 

Estes coeficientes são os que foram atribuídos 
no processo de simulação cujos resultados consti- 
tuem o presente trabalho. 

Eles foram escolhidos com base em razões téc- 
nicas que se imaginaram para o projecto em deter- 
minada fase do seu estudo. 

Obviamente, esses coeficientes podem ser alte- 
rados quer para experimentação do processo de si- 
mulação, quer para, no caso real, corresponder a 
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FIG. 1.5.1 - DISTRIBUIÇÕES TRIANGULARES DOS COEFICIENTES 


uma situação «de facto» referente ao conhecimento 
num determinado momento. 

Na nossa realização do processo de simulação 
vamos por certo alterá-los, 

Uma nota final sobre este aspecto das variáveis 
aleatórias do modelo aqui apresentado deve referir 
que nele não se entra em linha de conta com a alea- 
toriedade das reservas, o que, de certo modo está 
em contradição com o que se diz na descrição do 
projecto (ANEXO I — 5." página). 

E não será só este aspecto, por certo, que pode 
ser criticável no modelo apresentado. 

De resto, de todo este processo de simulação, o 
modelo matemático dificilmente traduzirá toda a rea- 
lidade, e será suficientemente genérico para abarcar 
todas as hipóteses, a menos que os custos e os tem- 
pos de realização não sejam considerados. 

Por outro lado, são estes aspectos (custo e 
tempo) de um modelo suficientemente genérico que 
levam à construção de um modelo com objectivos 
muito específicos, e cujo aproveitamento para objec- 
tivos diversos ou outro projecto, é geralmente im- 
possível. 


5.2 — Correlações entre variáveis 


Da observação dos esquemas representativos das 
distribuições triangulares atribuídas às variáveis, 
deve ressaltar a 2) à qual o valor mais provável 
atribuído não é representado pelo coeficiente 1, mas 
sim .79. 
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Nota-se ainda que o coeficiente aleatoriamente 
definido tem o valor máximo de 1 e por outro lado, 
a correcção que esse coeficiente irá provocar, vai 
incidir sobre uma variável já afectada de outro coe- 
ficiente (a variável 1). 


Trata-se na realidade dos valores de venda dos 
dois tipos de minério que o projecto prevê produzir. 
E tratando-se de dois tipos do mesmo minério, os 
valores de venda não são variáveis independentes 
entre si, mas também, uma não é função da outra. 


Isto é, considera-se absurdo que ao simular uma 
realização do projecto, ocorram para estas duas va- 
riáveis situações de valores absolutos incompatíveis 
com a análise histórica dos valores das mesmas va- 
riáveis. 


Trata-se portanto de uma estimativa de corre- 
lação que é introduzida no modelo. 


6 — PRÁTICA DO PROCESSO DE SIMULAÇÃO 


6.1 — Introdução de dados 


Para que o computador possa realizar a tarefa 
que através do programa ou programas lhe é indi- 
cada, torna-se necessário como todos sabem, forne- 
cer-lhe os dados (data) do problema a resolver. 


Mas também é sabido que essa data, tem uma 
certa ordem, tem um certo significado, e além disso 
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SSB REM LIMITES+*S20R A FORMA DE COEFICIENTES FARA 
PP? REM FRECO PERRLE 

1000 DATA +83s1+1,25 

10797 REM FRECO  SMALLS 

1100 DATA 165» Pol 

119% REM INVESTIMENTOS EM EQUIPAMENTOS E CAP. OPER,. FIXO (SOBRESSAL) 
1200 DATA .PSslod.3 

1299? REM INVESTIMENTOS EM CONSTRUCAO CIVIL 

1360 DATA «Psist.23 

1399 REM CUSTOS OFERACIONAIS E CAP. OPER, MOVEL 
1400 DATA +«Prlsisa 

1500 CHAIM SDCF2, 

2000 ENN 


READY 
SDCF'2. 


140 REM VALORES. DE BASE = ans ve ca 
149? REM PRECO DE VENDA DO PERRLE 48/T0N.CONC GFUBr) MARTA 
150 Pi=234 aa 

15% REM CUSTO OFER,. DIRECTO  ($8/TON, CONGO ado Td 

160 Ci=35,67 | a 

Lã? REM CUSTO DE TRANSPORTE FERROVIÁRIO (6/TON AM) 

179 REM NISTANCIA VIA FERREO € RM) 

180 K1=440 A 

189 REM CUSTO OPERACIONAL INHIRECTO (8/T0N.CONC) 

170 Ci=18.3 


4999 REM PRODUCÃO DE FPERELE 

000 DATA Lo DD. SDS DS PSP LP lvl Prios 
3099 REM PRODUCÃO DE SMALLS 
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v19PP REM GEOLOGIA E PROSPECAOQ 

5200 DATÁ 24105 15710,60312:2919040 

929? REM ESTULOS COORDENACÃO E OFTIMIZACAO 

5300 DATA 1944104 15710710+45+4951Lr1+1 

Bar? REM EQUIPAMENTO MINEIRO 

400 DATA 92,75,587,7550 

5499? REM EQUIPAMENTO DE TRANSE, Td, 

58500 DATA OrásB4 

559? REM ERIPAMENTO GERAL, 

3400 DATA ,857,2+0 

5479 REM EQUIPAMENTO DE TRATAMENTO E PARQUE 

5700 DATA 17,59157.5,0 

3799 REM EQUIPAMENTO DFICINAL 

3800 DATA 14759 15:75,0 | 

5099 REM EQUIFAMENTO DE DISTRIBUICAO DE ENERGIA 

3900 DATA 2.2+19.8B,0 

3999 REM  RARRAGEM E CENTRAL DE BOMBAGEM 

000 DATA 304740 

6099? REM RARRADEM LIFE LAMAS 

100 DATA 3,250 

819PP REM EXNIFICÇIOS INDUSTRIAIS 

6200 DATA 2072050 

629 REM. MOVIMENTAÇÃO DE TERRAS E ESTRADAS 

6300 DATA 251040 

4399 REM HABITACUES E EDIFICIOS SOCIAIS 

6400 DATA 30s1ãr5 

499 REM INSTALAÇÃO E TRANSFERENCIA DO TRANSP. -DE T.l, 
4500 NATA 0,0s520 ” 
6399 REM FERRO MADEIRA COMBUSTIVEIS LUBRIFICANTES FERRDSILICIO FLOCUL « 
400 DATA Br512 

499 REM CAPITAL OPERACIONAL MOVEL 

200 DATA 35+/100-30 

190000 END 
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tem, no caso geral, que ser introduzida em certos 
locais (números de linha) dos programas. 

No conjunto de programas (package) que cons- 
titui este processo de simulação apenas o primeiro 
e segundo programas necessitam de data a intro- 
duzir pelo interessado. Quanto ao terceiro programa, 
a sua data não é mais que o «sout-put: do segundo 
programa. 

Apresenta-se a seguir a data a fornecer, em 
ambos os programas e ela é suficientemetne expli- 
cita porque está precedida, para cada caso, da ins- 
trução REM. 


QUANTOS RUN ? 50 


ct b Y 


Se houver o cuidado de cotejar os valores da 
data do segundo programa, veremos que eles são os 
valores do «cash-flow básicos (Anexo 1), 

Quanto à data do primeiro programa, julgamos 
suficientemente explícito o conteúdo de 5.1. 


6.2 — Número de realizações 
Este aspecto encontra-se sumariamente abor- 


dado neste trabalho e é precisamente um dos as- 
pectos a desenvolver no futuro. 
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Creio que poderá dizer-se «quantas mais me- 
lhor», mas se neste caso o tempo não conta, (por 
cerca de cada cem «runs» um computador sofis- 
ticado pode demorar pouco mais de um minuto de 
unidade central), o preço já não será factor a des- 
prezar. 

O facto é que a experiência deve ser feita para 
cada caso, embora se saiba de antemão que os limi- 
tes superior e inferior, nos casos normais, têm valo- 
res de 600 e 300 respectivamente. 

As nossas experiências situaram-se quase sempre 
no limite inferior, e infelizmente a sobreposição de 
gráficos (que não são rigorosos uma vez que não 
são feitos em máquinas de desenhar) não é muito 
fácil, dadas as diferentes escalas de representação 
da variável. 

Uma primeira comparação foi feita com a inter- 
secção a 50 % de probabilidade, mas os seus resul- 
tados não são conclusivos. 


6.3 — Resultados 


Tendo adoptado trezentas realizações como va- 
lor do número de «runs» apresentam-se os gráficos 
dos dois outros indicadores escolhidos e calculados 
pelo modelo. 

Finalmente, só para exemplificar, fizeram-se va- 
riações nos coeficientes que traduzem a «insegurança 
do conhecimento» sobre duas variáveis. 
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Fig. 7.6 


O seu resultado (um só indicador — TDI) é ex- 
presso nos gráficos seguintes e o efeito dessa varia- 
ção analisado, mais uma vez e exclusivamente, 
através da intersecção a 50 % de probabilidade. 
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ANEXO | 


ANÁLISE ECONÔMICA 
DE EMPREENDIMENTOS MINEIROS * 


(Aplicação das técnicas do «Discounted Cash-Flown) 


1 — INTRODUÇÃO 


Os problemas tecnológicos cada vez mais com- 
plexos associados à legislação industrial dos nossos 
dias (por vezes não menos complexos que aqueles) 
e à tendência, *ou melhor, à necessidade de planear 
ou seleccionar projectos a longo prazo, originaram 
o desenvolvimento das técnicas para os estudos eco- 
nómicos. 

Por isso, e também pela consciencialização de 
que «lucro» só pode ter significado quando relacio- 
nado com o tempo e com o volume dos investimentos, 
a análise económica enveredou nos últimos anos pelos 
«métodos de desconto». 

Convenientemente elaborados, estes métodos for- 
necem os elementos de carácter económico necessá- 
rios às decisões, e posteriormente, em cada momento, 
constituem um aferidor da economia dos empreendi- 
mentos, o que, sem dúvida, é uma eficiente base de 
administração. 

A incerteza dos acontecimentos e portanto o 
«risco dos investimentos» não é contudo eliminado, 
mas a lógica dos métodos analíticos de desconto su- 
porta largamente a «tomada de decisão». 


Como em quaisquer outros métodos analíticos, a 
validade dos resultados depende da validade da infor- 
mação usada, mas até mesmo este aspecto pode ser 
investigado, isto é, podem ser determinadas as con- 
sequências finais (sensibilidades) das informações 
menos exactas. Resulta daqui, pois, uma definição de 
graus de sensibilidade ou, por outras palavras, a 
determinação das variáveis elementares que designa- 
remos por criticas, e, mais ainda, a possibilidade de 
estabelecer limites para os campos de variação das 
mesmas, 


Dentre os «métodos de desconto» usados em Aná- 
lise Económica, o «Discounted Cash-Flow» mostra-se 
eficiente na grande maioria das aplicações, sendo 
raros os casos em que as suas conclusões sejam me- 
nos significativas. 


2 — PRINCÍPIOS E TÉCNICAS DO «DISCOUNTED 
CASH-FLOW» 


2.1 — O Conceito de «Valor Actual» 


O poder potencial de crescimento do dinheiro ao 
longo do tempo, quando a ele se associa o conceito 
de juro, facilmente nos leva a concluir que as con- 


sequências das operações monetárias não são inde- 
pendentes do tempo. 


O conceito do «valor actual: de determinada 
quantia, refere essa quantia a uma data específica 
(na sua forma mais simples refere-se ao momento 
presente) e não é mais que a transformada daquela 
quantia num novo referencial de tempo. Note-se po- 
rém que o conceito do «valor actual» se baseia ex- 
clusivamente na razão atrás apontada, e nada tem 
que ver com a variação do poder de compra do di- 
nheiro, que normalmente é decrescente, 


A razão da transformação é da forma 


(1+r)º 
onde: 
n — são unidades de tempo. 


r— é o juro (expresso em forma decimal) 
por unidade de tempo. 


Da análise desta função verifica-se que o con- 
ceito do «valor actual», como já vimos, intimamente 
relacionado com o tempo, não pode abstrair-se do 
conceito de juro, a menos que o tempo seja nulo, o 
que origina a coincidência do «valor actual» com 
o valor efectivo. 


O «valor actual» (PV) calcula-se a partir do 
valor efectivo (V) dividindo este pelo factor de juro 
composto da unidade monetária (para o juro aim- 
ples por unidade de tempo, e unidades de tempo 
envolvidas na transformação). 


PV=V/(1+r)” 


Esta noção conduz efectivamente à confirmação 
de que o «valor actual» é igual ao valor efectivo 
quando o intervalo de tempo é nulo. 


PV=V/(1+r)º=V 


A aplicação mais corrente do conceito que su- 
mariamente tentamos expor, é a comparação de ver- 
bas cuja realização se processe em épocas diferentes. 
Tal comparação só se torna efectiva quando todas se 
encontram no mesmo referencial de tempo. 


Nesta aplicação, o juro utilizado para a trans- 
formação pode tomar qualquer valor diferente de 
zero, mas terá que ser o mesmo na transformação 
de todas as verbas que se pretendam comparar. 


Outra aplicação do conceito do «valor actual» 
é aquela que resulta da particularização deste mesmo 
valor, permitindo assim achar o juro a aplicar ao 
valor efectivo da verba para que esta se transforme 
em dado intervalo de tempo, no «valor actual» pre- 
viamente fixado. 


Resulta assim um novo conceito, como a seguir 
se expõe. 


* 1I Congresso Hispano-Luso-Americano de Geologia Económica. Madrid-Lisboa, 1971. 
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2.2 — Taxa de Desconto Interna 


O aspecto financeiro de um empreendimento ca- 
racteriza-se sempre por movimentos monetários em 
sentidos opostos, isto é, despesas e receitas, 

Quando estes movimentos monetários se agru- 
pam nos intervalos de tempo em que se processaram, 
a soma algébrica daqueles conduz a um valor final 
(saldo) próprio de cada intervalo de tempo, e que 
tanto pode ser positivo (receita) como negativo 
(despesa). 

Quer dizer, em última instância, um empreendi- 
mento normal terá nos intervalos de tempo da sua 
vida, saldos positivos e negativos. 

A apreciação do conjunto dos valores efectivos 
dos saldos só se torna realmente válida se, como 
vimos atrás, os situarmos no mesmo referencial de 
tempo, isto é, se calcularmos o «valor actual» de 
cada um deles. 

Chama-se «taxa de desconto interna» aquele 
juro que, aplicado aos saldos positivos, lhes confere 
«valores actuais» cuja soma é igual mas de sinal 
contrário à soma dos «valores actuais» dos saldos 
negativos quando descontados ao mesmo juro, 

Chama-se portanto «taxa de desconto interna» 
aquele juro que confere à soma algébrica de todos 
os «valores actuais», um valor nulo. 


2.3 — Recuperação do Capital Investido 


O tempo necessário para recuperar os investi- 
mentos feitos ou a fazer com determinado empreen- 
dimento, sem qualquer juro, foi sempre e ainda hoje 
é um indicador de certo significado. 

A análise económica através dos métodos de 
desconto, mais concretamente o «Discounted Cash- 
-Flow» na sua forma cumulativa não descontada, 
fornece automaticamente aquele tempo. Porém, o seu 
significado é diminuto, a menos que o empreendi- 
mento seja sujeito a alto grau de «risco». 

Este indicador, de facto, é independente dos 
acontecimentos posteriores à data que ele define, e 
como tal nada nos diz sobre a restante parte da vida 
prevista para o empreendimento, 

Já bastante mais significativos se tornam os 
tempos de recuperação dos investimentos quando cal- 
culados (actualizados) a uma taxa particular, por 
exemplo a taxa mínima a que a empresa espera ver 
recompensado o seu capital próprio. 

Este raciocínio leva-nos, aliás, ao conceito já 
definido da «taxa de desconto», que afinal não é 
mais que o juro a que se terá que descontar o «cash- 
-flow» para que o tempo de recuperação do capital 
investido seja igual ao tempo de vida do empreendi- 
mento. 


2.4 — Estruturação do D. C. F. 
As técnicas e os princípios gerais envolvidos no 


«Discounted Cash-Flow», sumariamente abordados, 
mostram uma relativa simplicidade. 
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A estruturação de um DCF requer porém a co- 
laboração de todos os sectores envolvidos no em- 
preendimento e só na medida em que estes se inte- 
gram nos objectivos, poderão fornecer o melhor do 
seu conhecimento e da sua sensatez. 

Só neste caso o modelo económico traduzido pelo 
DCF terá alta probabilidade de corresponder à futura 
realidade, e portanto tornar-se útil, na medida em 
que previu os acontecimentos e ou constitui um afe- 
ridor para eles. 

Vejamos pois, sumária mas especificamente, os 
mais importantes aspectos. 


2.4.1 — Administração 


Sob o aspecto da estruturação, a Administração 
precisa definir claramente os objectivos, quer no as- 
pecto qualitativo, quer ainda, se possível, no aspecto 
quantitativo. 

A orientação geral inicialmente definida não 
deve porém ser considerada imutável, mas constituirá 
a base dos pareceres prévios dos sectores técnicos, 
legais, fiscais e financeiros, a submeter posterior- 
mente à mesma Administração. 

O exame conjunto dos referidos pareceres poderá 
alterar ligeira ou profundamente a orientação inicial. 
Quer dizer, as diversas alternativas que, perante um 
objectivo, se pôem em cada um dos campos necessá- 
rios para o seu alcance, e que originalmente não 
eram por certo do conhecimento da Administração, 
podem alterar a orientação básica definida em pri- 
meira instância, e podem até alterar os objectivos 
finais. 

Também as incertezas de cada operação unitá- 
ria, mesmo nesta fase prévia, devem ser definidas 
no parecer de cada sector, e à Administração com- 
petirá apreciá-las no conjunto. 

Nota-se no entanto que deduções mentais ou 
mais apressadas acerca de situações hipotéticas não 
podem ser consideradas sem que o estudo económico 
seja totalmente recalculado. A interacção prove- 
niente de variações básicas é praticamente imprevi- 
sivel, e o seu estudo é tema de capítulo posterior, 
fora da «estruturação» que estamos a abordar. 


2.4.2 — Técnicas 


O grau em que as técnicas se encontram envol- 
vidas na estruturação de um «cash-flow» é variável 
de empreendimento para empreendimento, mas duma 
maneira geral é muito importante. 

O conhecimento das técnicas é susceptível de 
aperfeiçoamento e, muitas vezes, esse conhecimento 
pouco diz sobre o caso específico que se pretende 
tratar. 

O caso concreto de determinadas técnicas apli- 
cadas a determinados produtos é pois um processo 
moroso que normalmente se processa por estágios. 


Na medida em que os estágios de conhecimento 
aplicado se aprofundam, as margens de erro dimi- 
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nuem, e por isso mesmo a informação técnica que 
enformará os estudos económicos deverá, em está- 
Eios sucessivos, conter valores cada vez mais aper- 
feiçoados e com eles serem definidas margens de 
erro cada vez menores. 

Quer dizer, durante um período mais ou menos 
longo, o estudo económico irá sendo aperfeiçoado, 
preferivelmente por estágios, fundamentalmente es- 
tágios de conhecimento técnico, através do aprofun- 
damento de determinados caminhos, do estudo de 
alternativas ou ainda da optimização de conjuntos. 

Por seu lado, o estudo económico de cada está- 
gio pode ir definindo concretamente os campos que 
merecem estudo mais acurado, isto é, os campos de 
maior sensibilidade do empreendimento, como vere- 
mos mais adiante. 


2.4.3 — Vendas 


Compreende-se facilmente, embora por vezes esta 
realidade fundamental seja esquecida ou menospreza- 
da, que as vendas determinam o interesse e, por- 
tanto, a necessidade de um projecto. 

O ajuste de qualidade entre o que se sabe e 
pode produzir e aquilo que se é capaz de vender, 
é o primeiro tema determinante do interesse do 
empreendimento. 

O segundo tema do sector vendas refere-se ao 
ajuste entre as quantidades vendáveis e produzí- 
veis, e o terceiro aos preços de venda, mas entre 
estes dois últimos a interacção pode ser funda- 
mental. 

Também as previsões fundamentadas têm ca- 
rácter muito importante neste campo, mas elas 
devem ser bem distintas das realidades existentes. 

O «risco» total do empreendimento é, em al- 
guns casos, constituído por larga parcela prove- 
niente do campo das vendas. 

Os técnicas de «marketing» e as campanhas 
de publicidade podem constituir elementos impor- 
tantes na concretização das vendas, mas nem to- 
dos os empreendimentos são altamente sensíveis a 
tais manobras. 

A prospecção de mercados desempenha um pa- 
pel primordial neste campo, e o carácter evolutivo 
do passado permite em maior ou menor grau pre- 
ver o futuro, 


244 — Impostos e Finanças 


O campo dos impostos é, de um modo geral, 
altamente significativo, não só pelo valor que os 
impostos podem tomar, como também, e duma ma- 
neira geral, pela possível complexidade do pro- 
blema. Os aspectos financeiros e legais estão tam- 
bém relacionados com estes, pelo que a análise 
conjunta se torna necessária, 

Constitui no entanto um bom método de tra- 
balho, a análise do empreendimento pelos próprios 
méritos, isto é, manter a distinção entre o em- 
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preendimento em si mesmo e os benefícios que 
certas manobras legais ou financeiras possam tra- 
zer no campo dos impostos. 

Um outro aspecto de importância no campo 
dos impostos é o da depreciação dos bens envol- 
vidos no empreendimento. Sendo a maioria dos im- 
postos calculados na base do lucro anual bruto 
menos as depreciações e/ou perdas acumuladas, é 
corrente a situação em que as técnicas mais habi- 
tuais de depreciação não sejam as mais convenien- 
tes para fins fiscais, desde que estes permitam 
certa maleabilidade. 

Põe-se então o problema de ensaiar, dentre as 
possíveis, qual a técnica de depreciação mais bené- 
fica ao empreendimento. 

O recurso ao financiamento a médio e longo 
prazo é política comum nos empreendimentos in- 
dustriais. A razão fundamental situa-se nas limi- 
tações do capital próprio, donde aliás provém a 
expressão atrás usada «recurso ao financiamento». 

Ora, de facto, a utilização de financiamentos 
tem efeitos mais vastos que o simples «tornar pos- 
sível o empreendimento». 

Porém, hã a considerar a limitação natural dos 
financiamentos, pela sua relação com o capital 
próprio, a menos que exista a possibilidade de um 
caval», 

O efeito dos financiamentos nos resultados do 
empreendimento, pode prever-se em termos de 
«sensibilidade», de que falaremos mais adiante. 

As situações em que o efeito dos financiamen- 
tos é benéfico, nulo ou prejudicial, são facilmente 
definíveis em termos gerais. 

Mas o que efectivamente interessa é a previsão 
de uma das três situações atrás enunciadas e, mais 
ainda, do seu grau de sensibilidade, 


2.5 — O Valor do Elemento Humano 


A potencialidade que advém do conhecimento 
e do interesse que os elementos humanos podem 
pôr ao serviço de uma empresa, é um meio de 
realização do empreendimento não traduzível no 
campo dos estudos económicos. Contudo, o seu 
efeito sobre os resultados é extremamente impor- 
tante. 

Por isso, uma verdadeira «equipes tem que 
ser associada ao empreendimento, no seu estudo bá- 
sico, na sua optimização, na sua realizacão e final- 
mente na sua exploração. 

Só o esforço integrado de todos os elementos 
da «equipe» — administradores, técnicos, economis- 
tas, advogados, etc. — associado ao conhecimento 
básico dos fins a atingir e também dos meios ou 
técnicas de análise económica, pode prever os ele- 
mentos estruturais da mesma análise, 

Desse esforço, é altamente provável que re- 
sulte um modelo econômico cujos desvios da rea- 
lidade sejam de menor importância e, por isso, os 
objectivos finais que deram origem à decisão de 
acção, sejam alcançados. 
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2.6— O D. C. F. como Auxiliar de Administração 


Como órgão máximo de uma empresa, a Admi- 
nistração deve estar absolutamente segura das téc- 
nicas, méritos e deficiências do D. C. F.. Se assim 
não acontecer, então ou não tira proveito daquilo 
que o D. €. F. lhe pode proporcionar, ou, aceitando 
incondicionalmente as suas conclusões, delega nos 
peritos económicos o poder de decisão, que é prer- 
rogativa própria dessa mesma Administração, 

Nesse aspecto a Administração baseia-se fun- 
damentalmente no que a seguir se expõe. 


2.6.1 — Sensibilidades do Investimento 


O resultado global do D. C, F. não é igual- 
mente dependente de todos os componentes que o 
euformam., 

Leste modo, o estudo da variação global por 
efeito de variação específica de um componente 
pormite estabelecer o grau de sensibilidade do em- 
preendimento em relação a esse componente, 

Esta técnica define portanto graus de sensi- 
bilidade para os vários componentes, do que resulta 
uma ordem para aprofundamento dos conhecimen- 
tos e, portanto, a exigência de valores médios mais 
precisos e «limites de confiança» mais apertados. 
O mesmo é dizer que os componentes mais sensi- 
veis são os que necessitam de ser mais bem do- 
minados, não só nos estudos, como, posteriormente, 
na execução. 

Os estudos de sensibilidade podem estender-se 
a todos os campos, isto é, a todos os sectores con- 
tribuintes, mas mesmo dentro destes, aspectos es- 
pecíficos podem ter interesse para determinação 
das causas últimas. 


2.6.2 — Consequências dos Atrasos ou Para- 
gens 


O aspecto que aqui se pretende focar não é 
mais que um caso específico da sensibilidade. 

Nota-se no entanto que, na generalidade, os 
empreendimentos são muito sensíveis a estes fenó- 
menos, razão pela qual se justifica o desenvolvi- 
mento que também recentemente têm tido as téc- 
nicas de PERT e C. P. M.. 

A aplicação destas técnicas associadas ao 
D. C. F, conduzem por vezes a decisões perfeita- 
mente justificadas que optimizam o todo de um 
empreendimento, mas que, analisadas por si, pare- 
cem perfeitamente absurdas. 


2.7 — Dinâmica do D, C. F. 


Já vimos que o D, C F. não é mais que um 
modelo económico do empreendimento, e também 
vimos que esse modelo é, na sua estruturação bá- 
sica, prévio aos acontecimentos. 


248 


Se admitirmos que esse modelo económico foi 
construído com bases técnicas, económicas, finan- 
ceiras e legais sólidas, o modelo constituirá um 
ideal do desenrolar dos acontecimentos, 

Evidentemente que o desenrolar ideal dos acon- 
tecimentos vai conduzir a resultados que estavam 
integrados nos objectivos finais da empresa, pois 
foi essa a razão da decisão de acção. 

Assim que os acontecimentos se começam de 
facto a desenrolar, o modelo econômico (cash-flow) 
deve ser comparado com eles ou corrigido naquela 
parte que os factos consumaram. 

A evolução mais provável de alguns dos acon- 
tecimentos reais começa cedo a definir-se e a acção 
para corrigir os afastamentos do modelo ideal ou 
manter a identidade com ele, pode e deve iniciar-se 
logo de seguida. 

É este aspecto vivo de comparação, correcção 
e nova previsão entre modelo e realidade que se 
pretende expressar pela «dinâmica do D. €C. F.s. 


2.8 — O Uso de Computadores no Cálculo de D.C. F. 


O cálculo da «taxa de desconto interna» de 
um «cash-flow» é operação demorada e fastidiosa, 
muito especialmente se ele for irregular na totali- 
dade ou em grande parte da vida do projecto. 

Manualmente, hã que recorrer a tabelas, que 
em função de uma taxa de desconto estimada e 
do tempo (unidade de tempo) nos indicam o «valor 
actual» de uma unidade monetária, 

Procurando-se a «taxa de desconto interna» 
haverá que refazer os cálculos com diferentes ta- 
xas de desconto, de modo que o «valor actual» 
tenha sinal algébrico diferente em duas delas. 

A interpolação linear entre as taxas de des- 
conto estimadas e os «valores actuais» respecti- 
vos, dá uma primeira aproximação do resultado, 
mas porque, quanto mais afastados são os valores 
menor é a precisão da interpolação, torna-se neces- 
sário, em muitos casos, entrar num método ite- 
rativo. 

A programação dos cálculos baseia-se obvia- 
mente na fórmula que deu origem às tabelas, e 
não nestas, isto é, baseia-se no conceito de «valor 
actual» que se traduz por: 


do 
PV = > EA? 
n (1+r)” 
A iteração necessária para aproximação da 
«taxa de desconto internay processa-se pela intro- 


dução de um valor r', correcção necessária da «taxa 
de desconto inicial», que se traduz por: 


(1+) > vt” 
n 


> nV (1+r)" 
n 
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3 — AS CARACTERÍSTICAS GERAIS DO PROJECTO 


3.1 — Generalidades 


Admitamos que numa região tropical, a cerca 
de 400 kms do mar, são conhecidos diversos jazigos 
superficiais de minério de ferro, dispersos em área 
relativamente restrita, pouco acidentada e mesmo 
pouco arborizada, atravessada por rio de caudal 
razoável, 

A zona é praticamente virgem de população, 
embora a distâncias da ordem dos 30 kms existam 
estradas, caminhos-de-ferro e até uma povoação de 
razoáveis dimensões. 

Desta descrição, embora extremamente sumá- 
ria, parece inferir-se que existe um certo número 
de condições favoráveis à exploração de tais jazi- 
gos, mas em contrapartida, considerando que só a 
exportação é viável, aparece-nos uma longa distância 
ao mar. 

O conhecimento actual das reservas, no tipo 
de minério hematílico eluvial que, em primeira 
fase, se pretende explorar, é de 54 e 18 milhões de 
toneladas, respectivamente Reservas Certas e Pro- 
váveis, 

A probabilidade de novas ocorrências do mes- 
mo tipo de minério, dentro da área restrita já refe- 
rida ou na sua periferia próxima, é muito pequena, 
embora a hipótese de jazigos idênticos mas mais 
afastados, não seja de desprezar. 

O conhecimento actual das possibilidades e 
custos de exploração e tratamento físico para pro- 
dução de minérios vendáveis está já definido e pode 
fixar-se em 65% de rendimento médio em peso 
(para concentrados de 62 - 64% Fe) e 35667 + 
15 %/ton de concentrado sobre vagão (para mine- 
ração e transporte ponderado em relação ao T. Q. 
e rendimento médio em peso do tratamento). 

Admitamos ainda que no campo dos transpor- 
tes ferroviários já dispomos de proposta da enti- 
dade ferroviária assegurando o transporte desde 
a mina até ao porto de mar, parqueamento e car- 
regamento de navios, à taxa de $25/ton km sendo 
a distância global de 440 kms. 

Por outro lado, uma entidade produtora de 
energia eléctrica propõe-se fornecer a potência ne- 
cessária às barras de um posto de transformação, 
na mina, a custo que não interessa aqui especificar, 
mas que foi considerado nos custos operacionais 
directos já indicados. 

Ainda por outro lado, existe já um conheci- 
mento relativo à construção civil, movimentação 
de terras para esta, e equipamentos de mineração, 
tratamento, parqueamento, oficinal e geral, que é 
necessário executar e adquirir. Não se trata, po- 
rém, e compreensivamente, de um conhecimento 
exacto, 

A completar este panorama, dispõe-se de uma 
noção bastante correcta dos custos operacionais 
indirectos, do capital operacional fixo e móvel, e 
ainda dos impostos a que a actividade mineira está 
sujeita. 
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E finalmente, mas de modo algum menos im- 
portante, há o conhecimento do mercado internacio- 
nal dos minérios de ferro, as suas cotações actuais 
e passadas, assim como as evoluções da oferta e 
da procura que as originaram. 

Pode dizer-se, portanto, haver conhecimento 
razoável da situação em que o empreendimento, 
a realizar-se, se irá integrar, assim como um co- 
nhecimento razoável das reservas e dos custos ope- 
racionais. Ficam pois, com certa insegurança, os 
valores dos investimentos. 


3.2 — Estruturação do «Cash-Flow» Básico 


Tendo em atenção não só a limitação necessá- 
ria deste trabalho, como também os objectivos do 
mesmo, a estruturação do «cash-flow» básico 
(aquele que pretende analisar o projecto pelos pró- 
prios méritos) é feita com o mínimo de notas ex- 
plicativas, embora objectivamente, no Quadro I. 

Da sua análise e do anteriormente exposto so- 
bre o conhecimento de reservas, exploração, trata- 
mento físico, mercados, etc. infere-se que não es- 
tamos situados no tempo zero, mas sim no decorrer 
do ano 3, mais concretamente, nos primeiros meses 
desse ano. 

No decorrer dos três anos anteriores foram 
processados trabalhos e diversos estudos, consti- 
tuindo etapas de conhecimento sobre o possível 
investimento, a que corresponderam obviamente, 
com elevado grau de risco, certas despesas. 

Estamos pois admitindo que se torma neces- 
sária a análise (em termos económicos) do pas- 
sado, a previsão do futuro, e a tomada de nova 
decisão: 


— abandonar, em pura perda, o trabalho 
feito; 

— aprofundar os conhecimentos 
campos? ); 

— iniciar a execução do projecto. 


(em que 


É uma contribuição válida para a tomada de 
qualquer destas decisões, que vamos procurar ex- 
trair destes estudos, através dos indicadores se- 
guintes: 


— Máximo capital investido— (MCI— IM $) 

— “Tempo de recuperação do capital, sem 
juros — (TRC — anos) 

— Taxa de desconto interna (saldos bru- 
tos) — (TDIB — %) 

— “Taxa de desconto interna (saldos lígui- 
dos) — (TDIL — %) 


Para o modelo econômico básico os valores en- 
contrados são os seguintes: 


MCI  — 732.606 
TRO  — 87 
TDIB — 17,82 
TDIL — 16,84 


Anexos: Quadro I, Gráfico I, Folha Cálculo I 
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QUADRO | — LEGENDA 


1.1 — Produção de minério de alto teor «pebble». 
1.2 -— Produção de minério de alto teor «smalls>. 
1  — Produção total (106 ton). 


2.1 — Valor do minério (pebble — 64 % Fe) 
234800/tomn. 


22 — Valor do minério (smalls — 62% Fe) 
185800/ton. 


2 — Valor total do minério F. O. B. (10º $). 
3.1.1 — Geologia e prospecção. 

3.1.2 — Estudos, coordenação e optimização. 
3.1 —'Total de 3.1.1 e 3.1.2 (103 $). 


3.2.1 — Equipamento mineiro — 97 500/(3x)/v. resi- 


dual = 15% a 4 anos. 


3.2.2 — Equipamento de transporte T. Q. — 60 000/ 
'v. residual = 20 % a 12 anos. 


3.2.3 — Equipamento geral — 8000/(2x)/v. resi- 
dual = 10 % a 6 anos. 


3.2.4 — Equipamento de tratamento e parque — 
175 000/v. residual = 20 % a 12 anos. 


2) — Equipamento oficinal — 17500/v. residual 
= 20% a 12 anos. 


3.2.6 — Equipamento de distribuição de energia — 
22 000/v. residual = 20 % a 12 anos. 


Sé — Total do equipamento montado no local 
(3.2.1 a 3.2.6 — 103 5). 


3.3.1 — Barragem e central de bombagem — 
37 000/v. residual = 40 % a 12 anos 


3.3.2 — Barragem de lamas — 5000/v. residual nulo. 


3.3.3 — Edifícios industriais — 90 000/v. residual 
= 40% a 12 anos. 


3.3.4 — Movimentação de terras e estradas — 30 000/ 
/v. residual = 40 % a 12 anos. 


3.3.9 — Habitações e difícios sociais — 50000/v. 
residual = 40 % a 12 anos. 
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3.3.6 — Inst. e transferência do transp. T. Q. — 
20 000/(3x)/v. residual — 25 %. 


3.3 —'Total de obras (3.3.1 a 3.3.6— 103 5). 


3.4.1 — Sobressalentes e peças de desgaste — 
33 000/(10 % do v. do equip. inicial). 


3.4.2 — Ferro, madeira, combustível, Jlubrif., ferro- 
silício, floculante — 15 000. 


3.4 —'Total de capital operacional fixo (3.4.1 a 
3.4.2 — 103 8). 


3.5 — Capital operacional móvel — 180000 
(10º 4). 


3 — Investimentos totais. 


41 -—Custo operacional de mineração e trata- 
mento, ponderado em relação ao transporte 
de T. Q. e ao rendimento em peso dos tra- 
tamentos (35867/ton de concentrado coloca- 
do sobre vagão). 


42 — Transporte ferroviário e carreg. de navios 
($25/ton km — 440 km). 


4 — Custo total do minério F. O. B, (103 $). 
5 — Saldos brutos (103 8) — (2-3 - 4). 
6 — Saldos brutos acumulados. 


71 — Amortização máxima (50% dos saldos 


brutos). 
7.2 — Saldos brutos menos amortizações. 


ts — Imposto (%) = 
amort.)/1,8 x cap. social 
== 750000 X 103 $). 


100 x (saldos brutos — 
(capital social 


7 — Valor dos impostos. 

8  — Saldos líquidos (103 $). 

9 — Saldos líquidos acumulados, 

() —Os valores entre parêntesis representam, 


dentro da rubrica em que se encontram, va- 


lores negativos, 
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4 — ESTUDOS DE SENSIBILIDADE 


Pelas razões já apontadas, apenas são aqui apre- 
sentadas algumas hipóteses de variação dos factores 
elementares, e os estudos de sensibilidade, só num 
caso, consideram a variação simultânea de diversos 
factores. 

Estes estudos deverão, porém, num caso real, in- 
Cluir todas as hipóteses lógicas e prováveis. 

E pois comum a necessidade de análise de algu- 
mas dezenas de hipóteses. 


4.1 — Sensibilidade à Existência de Reservas Pro- 
váveis 


Se admitirmos que das reservas prováve;s dos 
jazigos apenas existe uma quantidade da ordem da 
quarta ou mesmo da quinta parte em que elas foram 
avaliadas, a vida do projecto reduz-se de dois anos. 
Como contrapartida, só se prevê que não haja neces- 
sidade de substituição da segunda frota de equipa- 
mento de exploração e transporte a pequena dis=- 
tância. 

Isto é, a produção total de minério vendável re- 
duz-se de 46,28 milhões de tons (9811 300 000800) 
para 37,88 milhões (8 051 500 000800) e o investi- 
mento total apenas sofre uma redução de 82 875 0005. 

O cálculo dos indicadores já citados formece os 
valores seguintes: 


MCT — 732.606 
TRC — 8,77 
TDIB  — 17,13 
TDIL  — 15,27 


Anexos: Gráfico II, Folha Cálculo II, 


4.2 — Sensibilidade a Investimentos Superiores 


Um dos aspectos citados anteriormente referia-se 
ao conhecimento pouco seguro do valor dos investi- 
mentos. 

Se admitirmos que os valores considerados no 
modelo básico para obras, equipamentos e sobressa- 
lentes, são ultrapassados em 15 %, o valor global dos 
investimentos ascenderá a 873 348 000800 contra 
774 825 000800 admitidos anteriormente. 

O cálculo dos mesmos indicadores fornece os va- 
lores seguintes: 


MCI — 817.446 
TRC — 9,33 
TDIB  — 15,33 
TDIL  — 14,17 


Anexos: Gráficos III, Folha Cálculo III. 


4.3 — Sensibilidade ao Financiamento 


Se admitirmos que o empreendimento tem possi- 
bilidades de ser financiado em 85 % dos nvestimea- 
tos originais em equipamento e sobressalentes, ao 
juro de 8 5% ao ano, amortizável em cinco anos (10, 
15, 20, 20, 20%) obteremos um financiamento no 
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valor de 55 300 000800, que vence um juro total de 
93 632 000500. 

Se além deste financiamento se conseguir um 
outro para 90 % do valor das obras originais, ao juro 
de 7% ao ano, também amortizável em cinco anos 
(10, 20, 20, 20, 20 %), teremos um financiamento no 
valor de 208 800 000800, que vence na totalidade, o 
juro de 50 813 000800. 

Em conjunto, os financiamentos ascendem por- 
tanto a 564 100 000800, vencendo juros no valor de 
144 445 000800. 

Acontece porém que o capital próprio investido 
até à data do hipotético financiamento apenas atinge 
181 000 000800, donde se deduz ser necessário um 
«aval», 

Calculando os indicadores já habituais, temos: 


MCI  — 276.262 
TCR  — 9,38 
TDIB — 23,80 
TDIL — 21,87 


Anexos: Gráfico IV, Folha Cálculo IV. 


4.4 — Sensibilidade ao Avanço da Produção 


Um outro aspecto a considerar, no caso de a 
tomada de decisão ser «iniciar a execução do pro- 
jecto», refere-se ao prazo que se interporá entre essa 
decisão e o início da produção, isto é, diz respeito à 
aceleração imprimida à elaboração do projecto defi- 
nitivo e respectiva execução. 

É a sensibilidade do projecto a um hipotético 
avanço na produção, que vamos apresentar, mas 
neste caso acrescentaremos a determinação do «va- 
lor actual» que tal avanço representa, uma vez que 
em valor absoluto o seu significado é praticamente 
nulo (visto que a produção está limitada às reser- 
vas, e se for iniciada mais cedo, acabará mais 
cedo). 

Através do valor actual do presente «cash-flow» 
descontado à «taxa interna de desconto» encontrada 
no modelo básico, podemos calcular prémios a atri- 
buir aos empreiteiros (obviamente menores que o re- 
ferido valor actual), ou, se for caso disso, comparar 
propostas de custos diferentes mas prazos de exe- 
cução também diferentes. 

No caso presente, constatamos que o máximo ca- 
pital investido é menor que o do modelo básico em 
63 462 000$00, que o tempo de recuperação do capital 
investido é também menor e, finalmente, que o valor 
actual do avanço da produção (cerca de 4 meses) 
atinge o valor de 31 647 000800. 

Os indicadores que temos vindo a tomar como 
referência, têm agora os seguintes valores: 


MCI — 669.144 
TRC — 8,45 
TDIB  — 19,22 
TDIL  — 17,90 
PV — 31.647 


Anexos: Gráfico V, Folha Cálculo V. 
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IO DATA -B6000,=1 4000, -40000,-27%01b60,-4742554451134175169,]El4la 
11 DATA 66604, |B0B25, | 63602, 157091, 162067,432039 

30 DATA D,6 

FLIN 


DISCTS 9:55 LISBOA D6/UT/TI] 


CASH FLOW,» NOT DISCOUNTEL: 
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4.5 — Sensibilidade Simultânea 


Para exemplo de sensibilidade simultânea enca- 
ramos três hipóteses (duas das quais já analisadas 
individualmente), todas elas agravando a situação do 
modelo básico: 


— Redução de Reservas Prováveis (4.1); 

— Investimentos Superiores (4.2). 

— Custos Operacionais Directos (Modelo 
Básico +15 %). 


Nestas circunstâncias, vamos encontrar para os 
indicadores habituais os seguintes valores: 


MCI  — 826.435 
TRC  — 1010 
TDIB — 11,93 
PO. — 108 


Anexos: Gráfico VI, Folha Cálculo VI. 


5 — CONCLUSÕES 
5.1 — Conclusões de Carácter Económico 


A fim de melhor analisar os indicadores obtidos 
para os diversos modelos, apresentam-se os mesmos 
sob a forma de Quadro, acrescidos do cálculo quan- 
titativo de sensibilidade quer de TDIL em relação a 
TDIB — S (L,/B)-—quer de TDIL em relação a TDIL 
do modelo básico — S (L/MB). 


Nestas circunstâncias, seriam possíveis as seguin- 
tes conclusões: 


i) investir capital num empreendimento que 
apresenta uma «taxa de desconto internas» 
aproximada de 16 % com a recuperação do 
mesmo a cerca de 6 anos (8,77 — 2,5) mos- 
tra-se razoavelmente aliciante. 


ii) Um problema fundamental põe-se no campo 
dos financiamentos porque eles são absolu- 
tamente necessários, e mais, pelo montante 

que atingem em face do capital próprio, re- 
querem forçosamente um «aval». 


iii) Uma acção esclarecedora dos aspectos ma- 
croeconómicos e sociais que o empreendi- 
mento porventura possa ter, constitui tema 
fundamental e imediato a desenvolver, no 
sentido de obter o referido «aval». 


iv) Os financiamentos, além de «tornarem pos- 
sível o empreendimento» trazem-lhe benefi- 
cios de carácter económico, aos quais é alta- 
mente sensível, e podem constituir a contra- 
partida de evoluções desvaroráveis em rela- 
ção ao modelo básico. 


v) Dentro dos campos analisados no presente 
trabalho, o empreendimento não se apre- 
senta com alto grau de «risco», pois que, se 
três condições desfavoráveis ocorrerem si- 
multâneamente (tendo em atenção que uma 
delas não é irreversível) a «taxa de desconto 
interna» ainda se situa na ordem dos 11 %. 


QUADRO Il 

Indicador 3.2* 41º 4,2% 4.3* 44% 4.5 

MES axa + 732606 7 32606 817.446 276.262 669.144 826435 
TRC 8,77 8,77 9,38 9,38 8,45 10,10 
TDIB 17,82 17,13 15,33 23,30 19,22 11,93 
TOIL cos 16,94 15,27 14,17 21,37 17,90 11,02 
S (L/B) —-8,30  —10,80 —7,50 —8,30 — 8,90 —T,80 
S(L/MB) .. Fã —6,55 — 13,28 +30,08 +9,55 —32,50 


* Referência do capítulo. 


As conclusões a extrair destes indicadores depen- 
dem fortemente da constituição, disponibilidades, ob- 
jectivos e panorâmica de outros investimentos da 
empresa que se propõe trabalhar o empreendimento. 

Procuraremos exemplificar com um caso, porven- 
tura o mais corrente, de uma empresa com disponi- 
bilidades reduzidas em relação ao empreendimento 
em causa, e portanto fundamentalmente preocupado 
com o projecto que acabamos de analisar. 
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vi) Tem significado importante o esclarecimento 
e o controlo dos investimentos, bem como 
a aceleração imprimida ao projecto defini- 
tivo e sua execução, para que o agravamento 
provocado pelo aumento daqueles e o bene- 
fício trazido pelo início de produção mais 
cedo, resultem em saldo nulo ou mesmo po- 
aitivo. 
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vii) O esclarecimento da existência total das Re- 
servas Prováveis não constitui uma preo- 
cupação dominante, do ponto de vista da em- 
presa, isto é, dentro do âmbito microeconó- 
mico, 


5,2 — Conclusões Gerais 


Os modelos económicos apresentados não são, 
de modo algum, todos aqueles a cuja análise seria 
conveniente proceder. 

Portanto, apenas tentâmos enquadrar alguns 
princípios de análise de investimentos num empreen- 
dimento mineiro, seguindo um caso hipotético e ana- 
lisando alguns aspectos, 

Recorda-se, por exemplo, que o preço de venda 
do minério foi sempre mantido constante, quer ao 
longo dos anos, quer ainda em todos os modelos; que 
os valores residuais das obras e equipamentos não 
se fizeram variar; que os custos operacionais toma- 
ram um valor médio não só em relação às operações 
de desmatação, limpeza de jazigos, métodos de ex- 
ploração e transporte de T. Q., como, e ainda, em 
relação ao rendimento em peso das operações de 
lavagem, concentração e calibragem, 

Obviamente não se desceu ao pormenor de con- 
siderar as variações citadas, quer pela índole do tra- 
balho, quer pela situação (em tempo) que atrás indi- 
camos, 

Em estágios superiores do conhecimento dos ja- 
gigos já seria conveniente pormenorizar, assim como 
separar as verbas dos saldos anuais (agora conside- 


radas «discretas:) em parcelas «discretas: e «<“con- 
tinuas>, 

Um caso concreto é pois suficientemente com- 
plexo mas fundamentalmente trabalhoso, para que a 
utilização de computadores (com preferência pelo 
sistema «Time - Sharing») traga proveitos inesti- 
máveis, 
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A energia nuclear nos países do Comecon 


(Continuação da pág. 208) 


A potência instalada nas centrais 
nucleares dos países do Comecon 
cifrava-se, em 1976, em 7500 MW, 
contra 1100 MW em 1971, Actual- 
mente, aproxima-se dos 10000 
MW, prevendo-se que em 1980 
atinja os 30000 MW, correspon- 
dendo a uma produção de 200 mi- 
lhões de MW,h. Assim, muito em 
breve, cerca de 10% do consumo 
de electricidade dos países do Co- 
mecon será de origem nuclear, per- 
centagem essa que, em 1990, atin- 
girã 40-66 qq na Bulgária, 26-35 % 
na Hungria, 13-15% na Polônia, 
22-26 %& na RDA, 10-14% na Ro- 
ménia, 30-42 % ma Checoslováquia 
e mais de 10% na URSS. 

A utilização da energia nuclear 
deu lugar à criação, para o sector, 
de dois organismos de carácter eco- 
nómico o Interatominstrument, li- 
gado aos aspectos de investigação 
e equipamentos, e o Interatomener- 
go para a organização da constru- 
ção das centrais. 

Hoje em dia, a maior parte das 
centrais construídas na Europa de 
Leste está equipada com reactores 
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VVER-d440 (440 MW) com sistema 
de arrefecimento por água; entre- 
tanto, prevê-se que muito em breve 
a sua capacidade atinja 1000 MW 
(VVER-1000), assinalando-se a 
construção de outros reactores a 
neutrões rápidos. 

URSS — A URSS possui cen- 
trais mucleares desde Junho de 
1854, data da entrada em serviço, 
em Obinsk, da primeira central 
nuclear do Mundo. Em 1975, a sua 
capacidade instalada atingia oz 4,7 
milhões de kW, ou seja 2,2% da 
capacidade total das centrais eléc- 
tricas do país, prevendo-se que em 
1980 a capacidde nuclear instalada 
se eleve a 13-15 milhões de kW, 
para um total de 67-70 milhões de 
KW. 

Das centrais nucleares projecta- 
das, duas serão edificadas em co- 
mum com os países europeus do 
Comecon que participarão no seu 
financiamento e construção de 
acordo com os planos estabelecidos 
em 1977-588 pela Comissão Perma- 
nente para a Energia do Comecon, 
Localizadas em Khmel'nickij e em 


Konstantinovka, na Ucrânia, terão 
uma capacidade individual de 4000 
MW. A Hungria, Checoslováquia, 
Polónia e URSS iniciaram já os 
trabalhos relativos à central de 
Khmel'nickij e da linha de alta 
tensão (750 kV) que a ligará à Po- 
lónia. 

Bulgária — A primeira central 
nuclear búlgara foi construída em 
Kozloduj com o auxílio da URSS 
e equipada com dois reactores com 
uma potência total de 880 MW. A 
sua produção diária de electricida- 
de é de 15 000-20 000 MWh, ou 
seja, cerca de 20 % de produção de 
electricidade do país. 

Encontram-se em construção 
dois novos reactores que terão 
igualmente uma capacidade de 
440 MW cada um, e que deverão 
entrar em funcionamento, respecti- 
vamente, em 1979 e em 1980, 
Nessa altura, a capacidade insta- 
lada totalizaráã 1860 MW e produ- 
zirá 50 % da electricidade do país, 

Hungria — A primeira central 
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RESUMO 


Formula-se na sua generalidade o problema da 
determinação do contorno final dum Desmonte a Céu 
Aberto. 

Descreve-se o algoritmo de Lerchs e Grossman 
com base em noções da teoria dos grafos. 

Introduz-se o método dos cones deslizantes apli- 
cado qdo algoritmo de Lerchs e Grossmam e dá-se 
um exemplo de aplicação. Descreve-se a estrutura 
básica de um programa que utiliza este método. 
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PROGRAMAÇÃO DE DESMONTES 


— Determinação do Contorno Final dum Desmonte a 
Céu Aberto («Ultimate Pit») 


1 — Introdução 


O problema que iremos abordar insere-se no 
quadro geral do planeamento de uma exploração a 
céu aberto. 

Em toda a sua generalidade o problema de de- 
terminação da corta optimizante é um problema de 
selecção a longo prazo, isto é, é um problema de cál- 
culo de reservas. 

Como tal faz intervir três tipos de estruturas, 
a saber: 
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FERNANDO HUMBERTO MUGE 
Universidade Técnica de Lisboa 


Instituto Superior Técnico 


ABSTRACT 


The problem of the determination of an ultimate 
pit contour in order to mazximize some objective 
fonction is formulated ím all its generality. 

It is described the Lerchs and Grossmann algho- 
ritm based in concepts of the graphs theory. 

It is introduced the moving cone method applied 
to the Lerchs and Grossmam alghoritm and tfs 
given an example of application. Tt is described the 
basico structure of a computer program which per- 
forms the above mentioned method, 


— A estrutura da mineralização (regionalização 
do espaço mineralizado); 

— a estrutura tecnológica que actua através dos 
processos de extracção e de concentração ou 
tratamento, sendo o seu impacto na definição 
da corta optimizante de natureza essencial- 
mente geométrica através da definição do 
tamanho, forma a implantação das unidades 
a extrair (bloco tecnológico de desmonte); 

— a conjuntura económica que implanta em cada 
época uma rede de mercados e de preços que 
condiciona fortemente as várias fases de um 
empreendimento mineiro. 


Quanto aos parâmetros económicos vamos supô- 
“los conhecidos através de um estudo prévio do mer- 
cado. 

Supunhamos ainda que ensaios laboratoriais nos 
permitiram definir o melhor método de exploração 
e/ou de tratamento do mineral ou dos minerais a 
explorar. 

Por outro lado é suposto que os trabalhos de 
reconhecimento executados nos permitiram atribuir 
a cada bloco tecnológico, com margem de erro acei- 
tável, estimadores das variáveis regionalizadas que 
intervém no cálculo de reservas (teores, rendimentos 
de processos de tratamento, etc...). Como se sabe a 
margem de erro será mínima se os estimadores refe- 
ridos forem obtidos por krigagem (v.g. Prática de 
Krigagem e Parametrização). 
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Posto o problema em toda a sua generalidade, 
o que se pretende é determinar o contorno final dum 
desmonte a céu aberto que conduza à maximização 
duma função objectivo. Um exemplo corrente duma 
tal função é a diferenca entre o valor total do minério 
e o custo total do desmonte do minério mais estéril. 


| B=AT m(x)-CT-C'T — para minério de teor 
! médio m(x) 


| = (ST — para estéril 


em que A: preço de venda do metal por tonelada 
de tal-qual 
C : custo do desmonte do minério por tone- 
lada de tal-qual 
C': custo do desmonte da tonelada de es- 
ril, 


Os modelos matemáticos que se descrevem a 
seguir visam alcançar o objectivo indicado tendo em 
atenção apenas restrições de tipo geométrico que, 
usualmente, são o máximo ângulo de talude admissÍ- 
vel e a profundidade as quais são determinadas pelo 
tipo de material a desmontar. 

Pretende-se pois determinar o valor de V que 
maximize : 


| B(x,w,z) dx dy dz 


* 1 


Onde V designa o volume associado ao contorno 
óptimo tal que em qualquer dos seus pontos a sua 
inclinação em relação a um plano horizontal não ex- 
ceda o ângulo a(x, y, Zz) fixado. (Angulo de talude 
máximo admissível). 


2 — Algoritmo de Lerchs e Grossmann 


O algoritmo de Lerchs e Grossmann permite, 
utilizando conceitos da Teoria dos grafos (v.g. Ane- 
xo) determinar o contorno final óptimo de uma ex- 
ploração a céu aberto, 


21 — Formulação do problema 


Considere-se o jazigo a explorar discretizado em 
blocos paralelipipédicos (blocos unitários) de volu- 
me V,. Note-se que a dimensão destes blocos está 
condicionada por aspectos tecnológicos (método de 
exploração adoptado) mas na maior parte dos casos 
a informação disponível não permite uma estimação 
suficientemente precisa de teor médio dos blocos so- 
bre os quais irá incidir a selecção futura [1]. 

O bloco a adoptar terá portanto que estar de 
acordo com a densidade de informação disponível 


(1) Repare-se que O par X,, x, é um par ordenado. 
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sob pena de o planeamento perder toda a validade 
visto basear-se em estimadores com margem de erro 
inaceitável. 

Posto isto pode-se associar a cada bloco de vo- 
lume V|/ um valor da função objectivo: B,. O que se 
pretende portanto é encontrar a combinatória de blo- 
cos que maximize: x B,, respeitando os constrangi- 


-— 


1 
mentos de ordem geométrica atrás referidos. 

Suponha-se que cada bloco i é representado pelo 
vértice x, dum grafo e considerem-se neste os arcos 
orientados (x,, x,) se o bloco j tiver pelo menos um 
ponto comum com o bloco i e se o desmonte de i 
implicar o de j. 

Obtém-se deste modo um grafo tridimensional 
orientado (1), 

G=(X, v) formado pelo conjunto X de vértices x, 
e pelo conjunto , de arcos (X, X,). 

Vê-se facilmente que qualquer contorno admissi- 
vel de desmonte, i.e. qualquer contorno que satis- 
faça as restrições de tipo geométrico impostas, é 
representado por um fecho do grafo G=(X, v). 

Raciocinando apenas para a fatia vertical repre- 
sentada na fig. que não implica perda de generali- 
dade, e notando que 


o grafo G define uma aplicação T' de X em X tal que 


segue-se, de acordo com a definição de fecho, que 


DP(x)=(x)=Y, i=1,2,..,5 
Dx) =(4) UP (x) 


== (X,) U (Xp X,» X,) a 


a "a 
Fe) =(g)urt)uUr()= 
(JU (Xp X, HJUMTOENUT E) UT x) = 
e (Xp Ho mia X,) E z, 
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Verifica-se portanto que qualquer fecho Y de massa total máxima, isto é, deseja-se conhecer o 


G =(X, y) é um contorno admissível. conjunto 

Associando a cada vértice x, de G=(X,,) O 
valor algébrico m, =v, — ec denominado por massa Y CIA tal que X 6Y => pr(tx) er 
do vértice x, pretende-se encontrar o fecho de G de e 


NORMALIZAÇÃO DE UMA ARVORE 
(EXEMPLO) 
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é máximo, 


2.2 — Algoritmo 


Inicialmente o grafo G indicado em 2.1 é amplia- 
do com um vértice fictício x e arcos (x, x,) de 


= E = a R o] 
modo a construir-se em G uma árvore inicial A mnor- 
malizada. 


O processo mais simples para construir A” con- 
siste em obter o grafo orientado (X, v ) em que y, 
é o conjunto de todos os arcos fictícios (x, x,). Tal 


= o 
grafo é a primeira árvore A. 


A árvore que é então transformada em sucessi- 


vas árvores A VA EDP A de acordo com determi- 
nadas regras até que não sejam possíveis novas 
transformações. O fecho máximo é então dado pela 
soma das massas dos vértices dum conjunto de de- 
terminados ramos bem identificados da árvore final 
(vidé Teorema em Anexo). 


Construída a árvore A entra-se num processo 
iterativo. A iteração de ordem k+1 transforma a ár- 


ALGORITMO DE LERCHS E GROSSMANN 


[EXEMPLO DE APLICAÇÃO) 


k k+1 
vor normalizada A na árvore normalizada A 4 


k 
Cada árvore A = (X, y) é caracterizada pelo seu 


k 
conjunto ,, de arcos e pelo conjunto Y dos seus 
vértices fortes. O processo iterativo termina quando 
Y é um fecho de G. 
A iteração de ordem k+1 contém os seguintes 
passos: 
k 
1. Construção da árvore normalizada A . 
2. Se existe um arco (x, X,) em G tal que 
k * 
x €Y ex € X-Y ir para 3, caso contrário 
passar a 7. 
3. Achar a raiz x do ramo forte que contém x. 
l 
4. Construir a árvore A substituindo o arco 


k 
(x, X,) de A pelo arco (X,, x). 


a 
k+1 
5. Normalizar A o que conduz a À 


6. Ir go número 1. 


“-k 
7. Fim. Y é o fecho máximo de G. 


Lerchs e Grossmann [4] além de justificarem 
os diferentes passos do algoritmo demonstram que 
se atinge o objectivo desejado com um número finito 
de iterações. 


ALGORITMO DE LERCHS E GROSSMANN 


[EXEMPLO DE APL ICAÇÃO) 
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Orgenigrama de blocos 


Leitura de dados 


Criação da árvore 
A, 


Pesquisa dos 


vértices x, é x 


| Na o | ns 
O 


Criação da 


arvore a 


Cara Estã Não 
= normalizada 


Normalizar 


Fig. 4 


3 — Método dos cones deslizantes aplicado ao algoritmo 
de Lerchs e Grossman (ALG.) 


cones deslizantes («moving cone methods») a vanta- 
gem de seguir os passos do ALG. 

Assim, o cone de um bloco (vértice) genérico 
X« é um sub-grafo fechado do grafo inicial G e o 
conjunto X,. dos seus blocos (vértices) é um fecho 
do grafo G. 

Por outro lado dada uma árvore A (sub-grafo 
de G) que contenha o cone do bloco x,., este identi- 
fica-se com o ramo R, podendo ser forte ou fraco 
cossoante a massa global (valor total) dos vértices 
que o constituem (cone forte ou fraco). O método 
identifica e conserva em memória os sucessivos cones 
fortes e fracos, 

Os sucessivos cones vão sendo transformados com 
adições e subtracções de blocos tal como na norma- 
lização das árvores construídas no ALG. 

Os sucessivos passos do método são apresentados 
na Fig. 8 e podem ser separados em quatro partes 
ou ciclos do programa. As quatro partes diferenciá- 
veis no método permitem a identificação de fechos, 
cones fracos e fortes e realizam as transformações 
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O primeiro grupo de instruções permite a identi- 
ficação dos blocos de minério existentes no depósito 
(blocos para os quais a função objectivo tem valor 
positivo). 

O segundo ciclo só intervém quando é encontrado 
um bloco de minério e pesquisa a existência de blocos 
suprajacentes que pertençam ao cone do bloco de 
minério não identificado no primeiro grupo de ins- 
truções e que não tenham sido previamente desmon- 
tados. Por definição, os blocos suprajacentes a um 
dado bloco de minério, são aqueles cujo desmonte é 
necessário para efectuar o desmonte desse bloco, 

O terceiro grupo de instruções efectua a norma- 
lização dos cones fortes e fracos. Este terceiro ciclo 
só intervém quando o bloco suprajacente é por sua 
vez membro dum cone fraco cujos vértices ou blocos 
constitutivos permanecem em memória. Neste caso, 
três situações podem ocorrer (vidé Fig. 5): 


1— O valor total do cone forte do bloco inicial soma- 
do com o do bloco suprajacente é negativo ou 
nulo. Neste caso, o cone forte não consegue 
suportar mais um bloco suprajacente e passa a 
fraco e o programa prossegue com a pesquisa 
de outro bloco de minério (1.º ciclo de instru- 
ções). 


2 — O valor total do cone forte do bloco inicial com 
o do bloco interseptado é positivo, mas o cone 
fraco a que pertence este bloco continua fraco 
mesmo sem ter que o suportar. Neste caso O 
bloco suprajacente é arrastado pelo cone do bloco 
inicial, desligando-se do cone fraco e o programa 
prossegue com o 2.º ciclo de instruções que pes- 
quisa novo bloco suprajacente, 

3 — O valor total do cone forte do bloco inicial com 
o do cone fraco a que pertence o bloco supra- 
jacente, é positivo. Neste caso ambos os cones 
podem suportar os blocos suprajacentes e os 
blocos do cone fraco passam a pertencer ao cone 
forte do bloco inicial. O programa prossegue com 
a execução do 2.º ciclo de instruções (pesquisa 
de blocos suprajacentes pertencentes ao cone e 
não desmontados). 

O quarto ciclo de instruções garante que O 
cone forte do bloco inicial seja um fecho do gra- 
fo inicial G. Isto porque os cones fracos intersep- 
tados pelo cone forte podem não ser completos, 
isto é, podem não atingir a superfície. No caso 
de tal acontecer o programa reentra no 2.º ciclo 
de instruções e a rotina continua através do 2.º, 
3.º e 4.º ciclos até que o cone forte passe a fraco 
com demasiados blocos suprajacentes ou até que 
o fecho esteja completo, Neste caso todos os blo- 
cos pertencentes ao cone forte do bloco inicial 
são desmontados. 

Em ambos os casos o programa prossegue 
com o 1.º ciclo de instruções até que todos os 
blocos de minério tenham sido examinados. 

Neste caso o programa executa ainda um 
5.º ciclo de instruções cuja finalidade é detectar 
a existência de cones fracos que, devido ao facto 
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de alguns dos seus blocos terem sido arrastados 
por cones fortes e desmontados, podem passar a 
fortes e serem eles próprios também desmon- 
tados. 


— Exemplo (v. g. Figs. 6 e 7) 


O 1.º ciclo de instruções inicia-se com a 
pesquisa de um bloco de minério. O bloco Té o 
primeiro a ser encontrado. Entra-se então no 2.º 
ciclo de instruções que permite procurar os blocos 
suprajacentes pertencentes ao cone do bloco 7 e 
ainda não desmontados, O primeiro bloco nestas 
condições é o bloco 3. Como o valor total da fun- 
ção objectivo após se ter acrescentado o valor 
do bloco 3 continua positivo, o 2.º ciclo de ins- 
truções prossegue a pesquisa de novo bloco supra- 
jacente e encontra o bloco 2. O cone passa a fraco 
com a adição do bloco 2 e regressa-se ao 1.º 
ciclo de instruções procurando novo bloco de mi- 
nério. 


O bloco de minério seguinte é 0 9. O 2.º ciclo 
de instruções adiciona os blocos 5 e 4 do cone, 
mas quando encontra o bloco 3 que pertence ao 
cone do bloco 9 e também ao cone fraco do blo- 
co 7, passa ao 3.º ciclo de instruções. O bloco 3 
é retirado do cone fraco e acrescentado ao cone 
forte do bloco 9 devido a este permanecer posi- 
tivo com o bloco 3 incluso, tal não acontecendo 
ao cone fraco do bloco 7 que permanece negativo 
sem o bloco 3 (vidé caso 2 fig. 5). O 2.º ciclo de 
instruções prossegue e como não há mais nenhum 
bloco pertencente ao cone do bloco 9, passa ao 
4,º cíclo de instruções. Como nenhum cone fraco 
foi acrescentado ao cone forte, a rotina afecta 
todos os blocos pertencentes ao cone do bloco 9, 
ao desmonte. 


O 1.º ciclo de instruções detecta o bloco 13 
e o 2.º ciclo adiciona o bloco & ao cone do bloco 13 
e seguidamente encontra o bloco 7 que pertence 
a um cone fraco (que inclui o bloco 2). O bloco 7 
é retirado do cone fraco e passa a pertencer ao 
cone forte do bloco 13 (3.º ciclo — 2.º caso). 
O 2.º ciclo prossegue até encontrar o bloco 2 que 
pertence a um cone fraco, O 3.º ciclo transforma 
o cone forte do bloco 13 em cone fraco ao adicio- 
nar-lhe o bloco 2 (1.º caso Fig. 5). 


O último bloco de minério é o 15. O 2.º ciclo 
prossegue até encontrar o bloco 7. Este bloco 
pertence ao cone do bloco 15 e também ao cone 
fraco que inclui os blocos 8 e 13. 


Neste caso o cone forte mais o cone fraco 
são positivos e o 3.º ciclo de instruções adiciona- 
-os (3.º caso v. £. Fig. 5). 
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Os blocos 6, 2 e 1 são adicicnados ao cone 
forte do bloco 13 à medida que o 2.º ciclo prosse- 
gue a pesquisa de blocos suprajacentes. Como 
houveram 2 cones fracos que interceptaram o 
cone forte do bloco 13, o 4.º ciclo de instruções 
vai pesquisar se algum deles está incompleto. 
Como ambos estão completos, todos os blocos 
referidos são afectos ao desmonte e o programa 
termina porque não se encontra em memória 
qualquer outro cone fraco que seja necessário 
recuperar (5.º ciclo de instruções). 


METODO DOS CONES DESLIZANTES APLICADO AO ALG 
[EXEMPLO DE APLICAÇÃO) 
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MÉTODO DOS CONES DESLIZANTES APLICADO AO ALGORITMO 
DE LERCHS E GROSSMANN 


(DIAGRAMA DE BLOCOS) 


INPUT: 
“MATRIZ DOS VALORES 


DA FUNÇÃO OBJECTIVO 
POR BLOCO 
- Nº TOTAL DE BLOCOS 


PESQUISA DE UM BLOCO DE FIM DO PROGRAMA 
MINERIO [VALOR POSI- - N OUTPUT: 


TIVO DA FUNÇÃO 08- EXISTE pá AaLG: il PLANTA DO DESMONTE 
JECTIVO| > BLOCO AÇO | ASSINALANDO OS NIVEI 
INICIAL ATINGIDOS 


PESQUISA DE UM BLO- PESQUISA DO 
CO SUPRAJACENTE QUE BLOCO-RAIZ DE- 
PERTE NÇA AO CONE DO SSE CONE FRA- 
BLOCO INICIAL E QUE | Co [NOVO BLO- 
NÃO TENHA SIDO DES- | CO INICIAL| 
MONTADO | 


TODOS 05 BLOCOS 
PERTENCENTES AO 
| CONE FORTE DO BLOCO 
INICIAL SÃO DES- 
MONTADOS 


O CONE FORTE DO BLOCO 
INICIAL FOI INTERSECTADO 
POR ALGUM CONE 
FRACO 


O CONE FORTE DO 
BLOCO INICIAL PASSA 
À CONE FRACO 


“ TODOS OS CONES FRACOS" 

QUE INTERSECTAM O CONE 

DO BLOCO INICIAL ESTÃO COM- 
PLETOS, ISTO E ATINGEM À 

| SUPERFICIE 


f O VALOR TOTAL DA FUN- 
CÃO OBJECTIVO PARA O, 
CONE DO BLOCO INICIAL É 
POSITIVO 


O VALOR TOTAL DO CONE 
DO BLOCO INICIAL COM O 
DO BLOCO SUPRAJACENTE, 


E POSITIVO 


O BLOCO SUPRAJA- 
CENTE PERTENCE A 
UM CONE FRACO 


FE ERP CORE" | /RENRANDO DSLOCO. 
| FORTE DOBLOCO INICIAL, | SUPRA JACENTE AO SEU 
SENDO OSEU VALOR CONE FRACO, ESTE PASSA 
| ACRESCENTADO DO VALOR A FORTE 
DO BLOCO SUPRAJACENTE 


O CONE DO BLOCO INICIAL ARRASTA 
TODO OCONE FRACO A QUE PERTENCE O 
BLOCO SUPRAJACENTE: 

-05 BLOCOS DO CONE FRACO PASSAM 
A PERTENCER AO CONE FORTE DO 
BLOCO INICIAL ; 

-D VALOR GLOBAL DO CONE FORTE E 

ACRESCENTADO DO VALOR DOS BLOCOS 
PERTENCENTES AO CONE FRACO 


O BLOCO SUPRAJACENTE 
E ARRASTADO PELO 
CONE DO BLOCO INICIAL | 
DESLIGANDO-SE DO 
CONE FRACO 


Fig. 7 
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ANEXO 


Vejamos alguns conceitos da teoria dos grafos 
que estão na base do algoritmo de Lerchs e Gros- 
smann: 

Grafo G de ordem n é definido como sendo a 
reunião de um conjunto numerável X de elementos 
X, (j=1, 2... n), chamados vértices de G, com um 
conjunto vy de pares de elementos de X, chamados 
arcos do grafo, Simboliza-se por G = (X, y). 

O grafo G define também uma aplicação multí- 
voca T de X em X tal que 


Pode simbolizar-se, portanto, também por G = 
= (X, T). 

Grafo parcial dum grafo G = (X, T) é um grafo 

'=(X, 4) em que para y X, 4X 67X, isto é, 

é um grafo com o mesmo número de vértices mas 
com deficiência de um arco, pelo menos, 

Sub-Grafo de um grafo G = (X, T) é um grafo 
G=(A, 7,) em que 


ACX A VXeA PT, X=(TXNA 


isto é, um sub-grafo de um grafo é um grafo parti- 
cular que contém apenas parte dos vértices e arcos 
do grafo total. 
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Caminho é toda a sucessão de arcos (u, us +.) 
tal que o extremo final de cada arco é o extremo 
inicial do seguinte, 

árvore é um grafo finito conexo sem ciclos e 
com pelo menos dois vértices. 

Árvore com raiz é uma árvore com um vértice 
particular — a raiz. Naturalmente que qualquer vér- 
tice pode ser raiz. 

Ramo e «ramificação» — O grafo que se obtéfh 
suprimindo um arco «|, numa árvore A com raiz tem 
duas componentes. A componente A, que não contém 
a raiz da árvore A é denominada ramo de A, 

A raiz do ramo A, é o vértice do ramo, adja- 
cente ao arco Hj 

Um ramo é uma árvore, e os ramos dum ramo 
designados por ramificações (twigs). 

Suponhamos que a cada vértice x, do grafo G = 
= (X, v) está associado o valor algébrico m,, deno- 
minado de massa do vértice x,. Admita-se que tal 
grafo é uma árvore A=(X, v). 

Cada aresta q, (arco «,) da árvore A define um 
ramo A, notado por A, = (X, v,). A aresta 4, de- 
signa-se por suporte do ramo A, e a massa M, de 
um ramo A, é a soma das massas de todos os seus 
vértices. A massa M, diz-se associada à aresta q, 
(arco u,) (nu, suporta a massa M,). 

Na árvore da fig. 2, 9-6-3-7 é um ramo o qual 
tem, p. ex., 6-9 como ramificação, 

Na árvore exemplificada o vértice x é o que se 
chama uma raiz fictícia pois podem referenciar-se, 
directa ou indirectamente, em relação a x todos os 
vértices da árvore, O vértice adjacente ao vértice 
fictício é tomado para raiz propriamente dito da ár- 
vore, (Vértice 1). 

Circuito é um caminho finito, formado pelos ar- 
COS uy, lo =, H DO qual o extremo inicial de q, coin- 
cide com o extremo final de q. 

Fecho de um grajo finito G = (X, T) é um fecho 


de T, isto é, é uma aplicação T: X -sX definida por 


A pt 
P(XD)=(A)UTGEDUT XD) U. 


De outro modo é um conjunto de vértices Y OC X 
tal que se 
po O T=-> PtX,) eY 


Prop. Se Y é um fecho de G = (X, T) então G = 
=(Y, T,) é um sub-grafo de G= (X, 7). 

Aresta dum grafo G=(X, v) é todo o conjunto 
de dois vértices X, e X, tais que 
(X, X,) € v e/ou (x, X)ev com X =X, 

Uma aresta é, portanto, um par de vértices liga- 
dos por um arco num sentido ou no outro, ou por 
dois arcos de sentidos opostos. 

Uma aresta é, geralmente, simbolizada por uma 
letra minúscula com um traço por cima, p. &, 


u=(X, X,) e o conjunto de todas as arestas por v. 
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Cadeia é toda a sucessão de arestas (yu, uy +) 
em que cada aresta u, está ligada a q, , por uma das 
suas extremidades e a u,, pela outra extremidade. 
Designa-se geralmente pelos vértices que contém. 

Ciclo é toda a cadeia finita que parte dum vér- 
tice X, e chega ao mesmo vértice. 

Ciclo é toda a cadeia finita que parte dum vér- 
tice X, e chega ao mesmo vértice. 

Grafo conexo é todo o grafo para o qual wy X,, 


X, (X, [X) existe uma cadeia de X. para X.. 


Arestas pe n 


Numa árvore com raiz fictícia x uma aresta u 
(ramo A,) é caracterizada pela orientação do arco q, 
em relação a x : 


— u, € uma aresta positiva (aresta —p) se o 
vértice final está mais afastado de x, que o 
vértice inicial. A, é então denominado de ramo 
positivo (ramo — p). 

—u, é uma aresta negativa (aresta —n) no 
caso contrário. Neste caso À, é um ramo ne- 
gativo (ramo —n). 


Classificação análoga se pode fazer para os ra- 
mos dum ramo. 


Arestas fe d 


Diz-se que uma aresta 4 (ou um ramo A,) é 
forte (f) se sendo positiva suporta uma massa estri- 
tamente positiva ou sendo negativa suporta uma 
massa nula ou negativa. No caso contrário a aresta 
(ou ramo) diz-se débil (d). 


vértices fe d 


Um vértice x, diz-se forte (f) se a cadeia entre 
x, e x, contém pelo menos uma aresta forte, caso 
contrário o vértice diz-se débil (d). 


Na fig. 1 são vértices fortes 2, 4, 8,3,6, 7 e 9. 


Árvore normalizada — Uma árvore diz-se norma- 
lizada se a raiz fictícia x é um vértice comum a 
todas as arestas fortes. 


Normalização duma árvore 


Qualquer árvore A dum grafo G pode ser nor- 
malizada substituindo os arcos (x,, x,) de arestas — p 
f por arcos fictícios (X,, X,) e os arcos (X,, X,) de 
arestas —n f por arcos fictícios (x, x,). Este pro- 
cesso é repetido até que todas as arestas —f te- 
nham x como um vértice comum a todas as suas 
extremidades, 


Lerchs e Grossmann [4] demonstram que: 
Prop. 1 Se um vértice x, pertence ao fecho má- 
vimo F (fecho de massa máxima) duma 
árvore normalizada A, então todos os 
vértices X, do ramo A. também perten- 
cem q y. 

O fecho máximo duma árvore normali- 
cada é o conjunto Y dos seus vértices 
fortes. 

Se num grafo orientado G” é possivel cons- 
truir uma árvore normalizada tal que o 
conjunto Y dos vértices fortes de A é um 
fecho de G', então Y é o fecho máximo 
de G'. (1) 


Teorema: 


(1) O grafo G' aqui considerado supõe-se ser o grafo Inicial G ampliado dum vértice fictício x, com uma massa nega- 
tiva e arcos fictícios (x,, x,) ligando x, a todos os vértices x,. Dado que x, não pode fazer parte de qualquer fecho máximo 
de G a introdução de arcos fictícios (x, x,) não afecta a resolução do problema formulado. 
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Programação de desmontes em lavra subterrânea 
(Estacionarização de parâmetros) 


1. BEJA COSTA (1) 
JORGE VALENTE (2) 


RESUMO 


4 exploração de minério de ferro da «Companhia Mineira do Lobito, 8. A. ER. L.» (Cassinga-Angola; 
6X 105 ton/ano de concentrados) foi programada, durante 6 anos, com sucesso na estacionarização de pará- 
metros, tais como o mimero de máquinas em operação nos desmontes, o teor dos concentrados e o respectivo 
rendimento em peso, assim como na repartição de finos (1/5 mm) nesses concentrados. Nesta programação 
de desmonte introduziram-se ainda alguns melhoramentos, tais como a optimização da recuperação em peso, 
com a introdução de um algoritmo computorizado e mantendo simultaneamente o teor em ferro dos con- 
centrados num certo valor pré-determinado; um outro melhoramento foi o estabelecimento de uma sequência 
de desmonte que define as combinações de blocos a serem desmontados simultaneamente ou sequencialmente, 
de forma a se assegurar as condições descritas anteriormente. 

Este método de programação é agora adaptado ao caso particular das minas subterrâneas e é demons- 
trada a forma de relacionar o cálculo dos diferentes tipos de equipamento destas minas pela primeira vez, 
segundo supomos. 

No final descreve-se um exemplo prático de aplicação do método; este exemplo refere-se à Mina do 
Moinho (Aljustrel) de «Pirites Alentejanas 8, A. R. La. 


SINOPSYS 


The exploitation of the iron ores of «Companhia Mineira do Lobito S. A. E. L> (Cassinga-Angola; 
bx 108 PPYF of concentrates) is based upon a sophisticated planning schedule which, for six wears, has 
been successful in keeping several pertinent parameter steady, such as the number of operating machines in 
the mine, the iron content and the weight recovery of the concentrates, às well as the proportion of fines 
(1/5 mm) im such concentrates. In this planning schedule some improvements was introduced, such as the 
optimization of the obtained weight recovery, with the introdution of a computerized optimization algorithm, 
and simultaneously keeping the iron content of the concentrates at a predetermined value; another impro- 
vement was the establishment of working criteria for the mine which define the combination of blocks to 
be simultaneously or sequentialhy treated in order to ensure the conditions described, 

Al this method is now, in the present text, adaptated to the particular case of the underground mines 
ond the 100y to relationate the calculation of the number of the diferent types of equipment used in that 
mines is shoum, we suppose, for the first time. 

In the end, a pratical example of the method is also described; this example concerns to the Mina 
do Moimho (Aljustrel) of the Company «Pirites Alentejanas, 8. A, R. Ls. 


1. OBJECTIVO-INTRODUÇÃO 


No I Congresso Hispano-Luso-Americano de Geologia Económica (I CHILAGE, Setembro 1971) J. Q. 


Rogado, J. M. Torres, J, €, Marques e D, T. Manso apresentaram duas comunicações complementares inti- 
tuladas: 


— «Optimização, por um método directo, da alocação de um minério a dois processos de tratamento»; e 
— «Um método de programação de desmontes». 


(1) — Eng.º de Minas (IST) — Director Técnico das Minas de S, João e Algares (Aljustrel), 
(2) — Eng.º de Minas (IST) — Engenheiro Chefe de Fundo das Minas do Lousal, 
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Na primeira, notando que a planificação da exploração de um jazigo toma, implicita ou explicita- 
mente, em consideração os resultados previstos no processamento do minério tal-qual, é anunciado e demons- 
trado o seguinte importante teorema: 

«Ao alocar-se a cada bloco um de dois tratamentos possíveis, a perda em metal é mínima quando se 
transferem para o tratamento que conduz ao concentrado mais rico, prioritariamente, os blocos para os 

a* Rt ZE, R, ZE, R,t, 
— , onde R=——— (rendimento em peso), t=-———— (teor) e à* é o opera- 


A*R 5 E, 3 E, R 


dor que traduz a transferência de um bloco genérico de um para outro método de concentração». 

Na segunda, demonstrado já o teorema anterior e exemplificada uma sua aplicação prática, os autores 
definem um método de programação de desmontes aplicável a jazigos regulares onde seja possível definir 
painéis ou blocos de exploração dos quais se conhecem o volume, o teor médio, a distância efectiva do respec- 
tivo centro de massas ao ponto de descarga, etc... O método consiste em estabelecer um critério simples e 
preciso de combinatória dos blocos que garanta, durante toda a exploração, a estacionaridade do centro de 
massa do jazigo e, portanto, a constância da frota de desmonte e transporte. O método permite ainda garan- 
tir a constância, dentro de certos limites impostos «a priori», de outro ou outros parâmetros pertinentes, 
como por exemplo o teor do concentrado obtido por tratamento do minério tal-qual daquele modo explorado. 


Este método foi utilizado nos últimos anos na programação de desmontes a céu aberto dos jazigos 
eluviais de ferro da área de Cassinga (Angola), da Companhia Mineira do Lobito SARL, embora adaptado 
às restrições locais, e nessa utilização e adaptação colaborámos, incluindo a sua total computorização e con- 
trolo e análise de resultados. Estamos, assim, plenamente conscientes das virtudes de tal método, a principal 
das quais é que não conduz à fixação rígida de um programa de exploração, visto que escolhida certa tra- 
jectória de desmonte esta pode ser frequentada arbitrariamente, desde que determinadas relações de massas, 
a combinar entre blocos, sejam respeitadas, 

Regressados de Angola e vindo ambos trabalhar em minas de lavra subterrânea, decidimos tentar 
adaptar os conhecimentos anteriores ao caso destas minas, generalizando o método, não esquecendo que é 
condição para sua aplicação que o domínio do jazigo onde se pretende usá-lo possa ser descretizado em 
blocos susceptíveis de serem desmontados mediante qualquer sequência temporal, isto é, a existência de cons- 
trangimentos sequenciais no desmonte dos blocos implica a sub-divisão do jazigo em domínios onde seja 
indiferente a ordem de exploração dos referidos blocos (domínios temporizados). 


quais é menor 


2. SIMBOLOGIA USADA 


Iremos usar, ao longo do presente trabalho, a seguinte simbologia: 


e, — caudal objectivo de tal-qual a extrair da Mina 

E, — tonelagem efectivamente extraída do piso i 

e=et+e tre, +... +e +... — repartição da tonelagem objectivo constante e, por cada piso da Mina 

r — rendimento real do poço de extracção 

E - rendimento horário do poço de extracção 

h — tempo de trabalho do poço de extracção 

e, — tonelagem objectivo de tal-qual a extrair do piso i 

A. Bi Co cw Mire — unidades de exploração do piso i 

la — tonelagem abatida em cada «corte» da unidade X, 

Ci — ciclo da pá auto-carregadora na unidade X, 

te. -— tempo de carga, descarga e manobras da pá auto-carregadora em cada ciclo 
L, “— distância de carregamento a considerar na unidade x, (ida e volta) 

v —- velocidade média de trânsito das pás auto-carregadoras 

P., -—- capacidade, em peso, de transporte na unidade de tempo da pá auto-carrega- 


dora na unidade X, 
- capacidade útil, em peso, da pá auto-carregadora 


rá 

y — rendimento global de utilização das pás auto-carregadoras 

n,, — número de pás auto-carregadoras trabalhando na unidade X, 

n, — número total de pás auto-carregadoras no piso i 

e, — tonelagem objectivo de tal-qual a extrair da unidade X, 

Z — constante que, por definição, é Z=LI e, lo. 

Ri, — reserva a explorar na unidade X,, em peso, de tal-qual 

m — número de actos de exploração para esgotar o piso 1 

A — distância média («centro de massas») ponderada de carregamento a consi- 


- 


derar para todas as unidades de exploração do piso 1 
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reserva a explorar no piso i, em peso, de tal-qual 


rendimento em peso médio ponderado do minério do piso i 


Hi 

c, teor de concentrados médio ponderado do minério do piso i 

Hix rendimento em peso do minério da unidade X, 

Ci, teor de concentrados do minério da unidade X, 

a ciclo de rolagem de uma composição carregada na unidade X, 

k, tempo de rolagem em carga da locomotiva no piso i transportando minério 

di da unidade X, 

t., tempo de rolagem em vazio da locomotiva no piso i em relação à unidade X, 

ti tempo total de manobras (nas toldas e na arreação) num ciclo de rolagem 

d.. distância de rolagem para o minério da unidade X, (ida) 

Va velocidade média de trânsito em carga da locomotiva 

v.=Kr, velocidade média de trânsito em vazio da locomotiva 

= constante que relaciona v, e v, 

pl, é capacidade, em peso, de transporte na unidade de tempo da locomotiva na 
unidade X, 

p capacidade útil, em peso, da locomotiva 

UR rendimento global de utilização das locomotivas 

N, número de locomotivas necessárias no piso i para se obter o objectivo e, 

Z constante que, por definição, é Z=E e, d, 

D, distância média («centro de massas») ponderada de rolagem a considerar 
para o piso i 

q capacidade útil, em peso, de cada vagona 

[8] - número de unidades de exploração em actividade de carregamento simulta- 
neamente no piso i 

N. número de ciclos da pá auto-carregadora necessários para completar uma 
composição de vagonas 

Es produção horária de uma composição de vagonas com minério da unidade X, 

N. número de composições de vagonas necessárias, numa hora, para a unidade X, 

N, número de vagonas em cada composição 

E, - ciclo do poço de extracção para o piso i 

H, profundidade do piso i 

V, velocidade de subida no poço de extracção 

Vi velocidade de descida no poço de extracção 

V=v,= velocidade considerada para o poço de extracção 

É tempo de carga e descarga no poço de extracção 

Po; capacidade, em peso, de extracção na unidade de tempo do poço de extracção 
em relação ao piso i 

Pp capacidade útil, em peso, do poço de extracção 


3. A UNIDADE DE EXPLORAÇÃO 


A estrutura do jazigo e a correlativa rede de acessos têm importância primordial na definição da 
unidade de exploração a considerar. Os autores dos dois trabalhos já citados usam os painéis ou blocos 
de exploração como sendo as suas unidades de exploração e a cada uma dessas unidades estão adstritos 
vários parâmetros, mas não impuseram constrangimentos temporais para a exploração daquelas unidades, 
pois consideraram que qualquer delas seria acessível para desmonte, em qualquer momento, independente- 
mente do avanço do desmonte nas unidades de exploração vizinhas. 


No desmonte subterrâneo usual de massas minerais segue-se actualmente o método de «câmaras e 
pilares» e, dentro de cada câmara, utiliza-se o «cut and fill». Atendendo precisamente à estrutura do jazigo, 
à sua rede interior de acessos e ao método de desmonte apontado, a unidade de exploração a considerar 
é o conjunto formado por duas câmaras contíguas (fig. 1). Poder-se-á perguntar qual a razão de uma 
unidade de exploração englobar duas câmaras e não apenas uma câmara. A resposta é fundamentada em: 

1.º) o maior investimento em equipamento verifica-se nas pás auto-carregadoras, pelo que há interesse 
em as manter devidamente ocupadas e evitar o seu sub-aproveitamento; para tal programa-se o desmonte 
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das duas câmaras da unidade de exploração de forma complementar, o que é possível, por exemplo, segundo 
o esquema: 


CAMARA A 


preparação 
perfuração 
ciclo comple- | disparo 
to de traba- | carregamento 


lho entulhamento 1 uso da pá 
entulhamento 
preparação 
perfuração 
(CAMARA A) (CAMARA B) disparo 
P né o o "” 
alxada para passagem de pessoal, 
ar, agua e materiais CAMARA B 
t tw s 
4 tolda (ponto de descarga do mi- carregamento 
nerio) ) uso da pá entulhamento 
entulhamento 
Ski ci é no . me reparação 
ES] chamine de entulhos e ventilaçao sipbi tera 
perfuração 
(existe so no "ceu" disparo ciclo completo 
a carregamento de trabalho 
“19. 1 
R uso da pá entulhamento 
Corte horizontal de uma unidade de exploração normal entulhamento 


(com 2 câmaras) 


2.º) Com o avanço do desmonte, a certa altura tem-se necessidade de proteger o equipamento exis- 
tente na câmara, o que é feito retirando-o para a outra câmara da unidade de exploração através da sua 
comunicação, aberta entre os dois pilares que dividem a unidade de exploração; todo esse equipamento fica- 
ria então inactivo, pelo que o respeito pelo apontado em 1.º tem a vantagem de ele estar protegido dos dispa- 
ros na câmara vizinha, para além de não ter equipamento de carga e transporte parado. 

No entanto não se deve esquecer as franjas do jazigo, que podem obrigar a existência de um número 
impar de câmaras no mesmo piso da Mina, Neste caso pode-se escolher a melhor de duas alternativas: 


1.º) a unidade de exploração ficará constituída 
zona esteril a muro por uma única câmara (fig. 2) mas há que realçar o 
aumento de custos devido à obrigação da existência 


de um maior número de toldas por área a desmontar 
e ao equipamento ter de ser protegido na zona aberta 
em estéril sem ficar em actividade; e 


2.º) a unidade de exploração ficará constituída 
por três câmaras (fig. 3), funcionando uma câmara 
como complemento de duas outras câmaras; nos casos 
de franja normais esta solução é preferível à da cons- 
tituição de uma unidade de exploração apenas com 
uma câmara. 


Fig. 2º 


Corte horizontal de uma unidade de exploração só com 
1 câmara 
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Fig. 3 


Corte horizontal de uma unidade de exploração com 3 câmaras 


Quanto ao problema da quantidade de equipamento tem de se notar que não é ele verdadeiramente 
que define a unidade de exploração, sendo adiante abordado o modo de determinação da quantidade de máqui- 
nas carregadoras e/ou transportadoras necessárias, apenas sendo certo que em cada unidade (quer tenha 1,2 
ou 3 câmaras) terá de haver o mínimo de uma pá auto-carregadora, podendo, no entanto, haver maior 
número de pás na frota de uma unidade de exploração. 


4. INSERÇÃO NO JAZIGO DAS UNIDADES DE EXPLORAÇÃO 


No caso geral um jazigo em exploração é dividido por pisos e cada piso em câmaras. Dada esta 
definição vamos considerar cada piso dividido nas nossas unidades de exploração, tendo a figuração esque- 
mática da fig. 4. 


SUPERFÍCIE 

Admite-se, como é usual, que se pretende 
extrair para a superfície caudal constante de miné- 
rio tal-qual explorado no jazigo (e, = constante), 
O poço de extracção terá, um rendimento real defi- 


nido por: 
Ex y Ez Ei 
r= + pe = ) 
e. e, e. e | 


por: 


No entanto, ao apontar-se uma repartição 

(e =e +86 + + e, + ...) da tonelagem objec- 
POÇO DE tivo constante e por cada piso da Mina, para se 
EXTRACÇÃO poder definir os rendimentos do poco de extrac- 
ção, é porque uma outra programação prévia já 
determinara, para o tempo h, qual o objectivo da 
extracção em cada piso da Mina. Então, vamos 


A, BE, €C,... — unidades de exploração, 
B, — unidade de exploração B do piso Y. 


Fig. 4 considerar o piso genérico i, do qual se pretende 

Esquema de um jazigo vulgar extrair, durante o tempo h, a tonelagem e,, tonela- 

gem esta a obter das diversas unidades de explo- 

ração do piso (A, B, C, ..., X, ...); cada uma destas unidades de exploração é considerada como uma uni- 


dade de trabalho (descrição justificada em 3), logo a cada unidade de exploração x. terão de ser adstritos 
os vários parâmetros e considerar e devidamente simbolizados pela colocação dos índices i. 
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